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ABSTRACT 

We assess the detectability of baryonic acoustic oscillations (BAO) in the power spec- 
trum of galaxies using ultra large volume N-body simulations of the hierarchical clus- 
tering of dark matter and semi-analytical modelling of galaxy formation. A step-by- 
step illustration is given of the various effects (nonlinear fluctuation growth, peculiar 
motions, nonlinear and scale dependent bias) which systematically change the form 
of the galaxy power spectrum on large scales from the simple prediction of linear 
perturbation theory. Using a new method to extract the scale of the oscillations, we 
nevertheless find that the BAO approach gives an unbiased estimate of the sound hori- 
zon scale. Sampling variance remains the dominant source of error despite the huge 
volume of our simulation box (= 2.41/i _3 Gpc 3 ). We use our results to forecast the 
accuracy with which forthcoming surveys will be able to measure the sound horizon 
scale, s, and, hence constrain the dark energy equation of state parameter, w (with 
simplifying assumptions and without marginalizing over the other cosmological pa- 
rameters). Pan-STARRS could potentially yield a measurement with an accuracy of 
As/s = 0.5 — 0.7% (corresponding to Aw « 2 — 3%), which is competitive with the 
proposed WFMOS survey (As/s = 1% Aw w 4%). Achieving Aw < 1% using BAO 
alone is beyond any currently commissioned project and will require an all-sky spec- 
troscopic survey, such as would be undertaken by the SPACE mission concept under 
proposal to ESA. 
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1 INTRODUCTION 



The discovery that the rate of expansion of the Uni- 
verse is apparently accelerating was on e of the key ad- 
vances in physical cosm ology in the 1990s (|Riess et al.l 19981 ; 
IPerlmutter et af]|l999l ) . Understanding the nature of the dy- 
namically dominant dark energy, which is believed to be re- 
sponsible for this behaviour, is one of the biggest challenges 
now facing cosmologists. 

Over the past decade our knowledge of the basic cosmo- 
logical parameters, which describe the content of the Uni- 
verse, its expansion history and ultimate fate has improved 
tremendously. This progress is the result of advances on two 
fronts: the advent of datasets which have provided fresh 
views of the Universe with unprecedented detail and the 
development of the theoretical machinery required to in- 
terpret these new measurements. Currently, the values of 
many cosmological parameters are known to an accuracy of 
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around 10% (albeit with caveats regarding degeneracies be- 
tween certain combinations of parameters and also regarding 
the precise number of parameters that are a llowed to vary 
in th e cosmological model; see, for example, ISanchez et al.l 
l2006t) . 

The cold dark matter (CDM) model has emerged as the 
most plausible description of our Universe. In the most suc- 
cessful version of this model, more than 70% of the density 
required to close the Universe is in the form of dark energy. 
Currently, there is no model which can reconcile the magni- 
tude of the dark energy component with the value expected 
from particle physics arguments. A simple phenomenological 
description of the dark energy is provided by the equation 
of state that relates its pressure, P, and density, p, which is 
encapsulated in the parameter w = P/pc 2 . If the dark en- 
ergy has the form of the cosmological constant, w = —1. The 
indications are that the dark energy now h as a form close to 
that expected for a cosmological constant (|Riess et al.ll200j : 
ISanchez et al.ll2006l ). However, in the absence of a theoreti- 
cal model for the dark energy, it is possible that the equation 
of state could depend on space and/or time. 

A whole range of experiments and surveys is being 
planned which number amongst their goals determining 



2 Angulo et al. 



the equation of state of the da rk energy as a func tion 
of redshift (for a discussio n, see lAlbrecht et al.l 120061 and 
IPeacock fc Schneiderl |2003 ). Several techniques are being 
considered, which are sensitive to the influence of the dark 
energy on various features of the cosmological world model. 
These include the Hubble diagram of Type IA supernovae, 
counts of clusters of galaxies, the weak gravitational lensing 
pattern of faint galaxies and the measurement of the bary- 
onic acoustic oscillation scale in the matter distribution as 
a function of redshift. The measurements and data analy- 
sis required to obtain useful constraints on the equation of 
state parameter are so demanding, and so open to poten- 
tial systematic errors, that it is necessary to pursue as many 
different avenues as possible. 

In this paper, we focus on the test using the baryonic 
acoustic oscillations (BAO). The BAO is the name given 
to a series of peaks and troughs on scales on the order of 
100 h^ 1 Mpc, imprinted on the power spectrum of matter 
fluctuations prior to the epoch of last scattering, when the 
matte r and radiation com ponents of the Universe were cou- 
pled (|Peebles fc Yul Il970f) . The BAO are the counterpart 
of the acoustic peaks seen in the power spectrum of the 
temperature of the cosmic microwave background radiation, 
though they have a different phase and a much smaller am- 
plitude dSunvaev fc Zeldovichlll970l;|Press fc Vishniadl 19801; 
Hu fc Sugivama 1996; Ei senstein fc Hulll99Sl ; iMeiksin etafl 



19991 ). The wavelength of the BAO is related to the size 
of the sound horizon at recombination. This does not de- 
pend on the amount or nature of the dark energy, but on 
the physical density of matter (Q m h 2 ) and baryons (Qf,h 2 ). 
Given the values of these parameters, for example, from the 
cosmic microwave background or large scale structure data, 
the sound horizon scale is known and can be treated as a 
standard ruler. The apparent size of this feature in the power 
spectrum of galaxies or galaxy clusters does depend on the 
dark energy and its equation of state th rough the angular di- 
ameter distance-redshift r elation (e.g. iBlake fc Glazebrookl 
l2003l ; lHu fc Haimanlliooj) 

BAO in the galaxy distribution were first glimpsed in 
the e arly stages of the " 2-degree-field galaxy redshift sur- 
vey" (|Percival et aL|[200ll ) and finally detected in the power 
spectrum of the completed 2dFGRS (|Cole et alJl2005T ). The 
equivalent feature, a spike, was also found in the corre- 
lation function measured from the luminous red galaxy 
(LRG) sample of the Sloan Digital Sky Survey (SDSS) 
I Eisenstein et al.1 120051 ). Cole et al used the BAO to con- 
strain the parameter combination (fijyi/fib, Qm) (where SIm 
and f2b denote the matter and baryon density parameters 
respectively). Eisenstein et al used the location of the spike 
in the correlation function to constrain the absolute distance 
to the median redshift of the SDSS LRG sample and hence 
constrained the value of Om- Hiitsi (2006a,b) carried out a 
power spectrum analysis of a similar LRG sample, and com- 
bined this measurement with other datasets to constrain 
the values of cosmological parameters. More recently, the 
BAO have been extracted from the power spectrum mea- 
sured from a much larger sample of SDSS LRGs to constrain 
n M and Oh/Orr, f iTegmark et all 120061; IBlake et al.1 [20071 ; 
iPadmanabhan et al.l 120071 ; iPercival et all 120071 ). To date, 
measurements of the BAO have only yielded constraints on 
the dark energy equation of state when combined with other 
datasets, such as the spectrum of temperature fluctuations 



in the microwave background or when restrictive priors have 
been adopted on certain parameters, such as the Hubble con- 
stant. 

The bulk of the work in the literature on the 
usefulness of the BAO has relied upon linear pertur- 
bation theory to assess the detectability of the fea- 
tures and to forecast the errors on the recovered value 



of w ( Blake fc Glazebrookl 120031; iHu fc Haimanl 



Glazebrook & Blake 2005; Blake & Bridle 2005; Blake et al. 



20061 ; iParkinson et al.ll2007l ). There are, however, a range 



of dynamical and statistical effects which can alter the 
appearance of the power spectrum relative to the lin- 
ear theory prediction, even on t he scale of the BAO , 
which we review i n this paper (ISeo fc Eisenstein! 120031; 
I Angulo etahl 120051; ISpringel et all 120051 ; ISeo fc Eisenstein! 
120051 ; Eisenstei n et al.l 20061 ). Some simulation work has 
been done to study these effects, mostly using computa- 
tiona l cubes of side 500 ft " 1 Mpc (ISeo fc Eisenstein! 12001 
l2005l ; ISpringel et ail 120051 ; lEisenstein et all 120061 ). These 
are only a small factor (2-3) bigger than the scale of the 
fluctuations of interest. Calculations with small boxes are 
subject to large sampling fluctuations and may even miss 
some features of the nonlinear growth of large scale fluc- 
tuations throu gh the absence of long wav elength density 
fluctuations (|Crocce fc Scoccimarrd l2006al ) . Very recently, 
larger simulation volumes have been used, of around a cu- 
bic gigaparscc and larger (ISchulz fc Wh ite 2006; Huff et al.l 



120061 ; lAngulo et all 120051 ; iKoehler et al] 120061 ). However, 



such studies have tended to have relatively poor mass resolu- 
tion, making it difficult to model galax ies without resor ting 
to simplified biasing prescriptions (e.g. ICole et al.iri998h . 

Given the significant commitment of resources required 
by the proposed galaxy surveys and the level of precision 
demanded by the BAO approach, it is imperative to en- 
sure that accurate theoretical predictions are available both 
to help in the design of the survey strategy and to extract 
the maximum amount of information from the observations. 
This is a tough challenge computationally, because it re- 
quires ultra-large volume N-body simulations with sufficient 
mass resolution to identify the haloes likely to host the 
galaxies to be seen in the surveys, and a realistic model 
to populate these haloes with galaxies. 

In this paper, we use a combination of suitable N-body 
simulations and a semi-analytical model of galaxy formation 
to assess the visibility of the BAO. In Section 2, we describe 
the suite of N-body simulations used and outline the semi- 
analytical model. Section 3 gives a blow-by-blow account of 
how the power spectrum changes relative to the simple pre- 
diction of linear perturbation theory, as additional layers of 
realism are added to the modelling, starting with dark mat- 
ter and ending with galaxies. We set out our approach for 
constraining the dark energy equation of state in Section 4, 
and present our results in Section 5. We give our conclusions 
in Section 6. 



2 METHOD 

In this section, we introduce the theoretical tools used to 
produce synthetic galaxy catalogues. First, we describe the 
N-body simulations (§2.1) which consist of a high resolu- 
tion run (§2.1.1) and an ensemble of lower resolution runs 
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Table 1. The values of some of the basic parameters used in 
the simulations. The columns are as follows: (1) The name of 
the simulation. (2) The number of particles. (3) The mass of a 
dark matter particle. (4) The softening parameter used in the 
gravitational force. In both cases, the length of the computational 
box is 1340 /i~ 1 Mpc, and the same cosmological parameters are 
used, as given in Section 2.1. 
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Figure 1. A test of the choice of starting redshift used in the N- 
body simulations. The upper panel compares the power spectrum 
measured at z = 15 in the BASICC when the simulation is started 
at z = 63 (dashed red curve) and at z = 127 (solid blue curve). 
The power spectra plotted in the upper panel have been divided 
by the linear perturbation theory prediction for the dark matter 
power spectrum at z = 15. The lower panel shows the ratio be- 
tween the power spectrum measured from the simulation started 
at redshift 63 to that measured from the run which started at 
redshift 127. 

(§2.1.2). Next, we discuss the measurement of power spectra 
from discrete distributions of objects and use the ensemble 
of low resolution simulations to estimate the errors on the 
power spectrum measurement (§2.1.3). In the second part 
of this section, we explain how a galaxy formation model 
is used to populate the high resolution N-body simulation 
with galaxies (§2.2). 

2.1 N-Body Simulations 

The N-body method is a long-established computational 
technique which is used to follow the growth of cosmological 
structures through gravitational instability (see, for exam- 
ple, the reviews by Bertschinger 1998 and Springel, Frenk & 



White 2006). Our goal in this paper is to simulate the forma- 
tion of structure within a sufficiently large volume to follow 
the growth of fluctuations accurately on the scale of the 
BAO, and with similar statistics for power spectrum mea- 
surements to those expected in forthcoming surveys. At the 
same time, we require a mass resolution which is adequate 
to identify the dark matter haloes likely to host the galaxies 
which will be seen in these surveys. To achieve these aims, 
we use a mem ory-efficient version of the GADGET-2 code of 
ISpringell (120051 ). which was kindly provided to us by Volker 
Springel and the Virgo Consortium. 

We use two types of calculation: a high resolution simu- 
lation, labelled the "Baryon Acoustic Simulation at the ICC 
or BASICC, which is able to track galactic haloes, and an en- 
semble of lower resolution simulations, labelled L-BASICC, 
which we use to study the statistics of power spectrum mea- 
surements on large scales. Here, we describe some of the 
common features of the simulations, before moving on to 
outline specific details in §2.1.1 and §2.1.2. 

We adopt a ACDM cosmology with th e same param- 
eters used in the Millennium Simulation (|Springel et al.l 
2005), which are broadly consistent with the latest con- 
straints from the cosmic microwave backg round data and 
large scale structu re measurements (|Sanchez et al.l 120061 : 
ISpergel et al.|[2007l '). The values of the parameters are: the 
matter density parameter, Qm = 0.25, the energy density 
parameter for the cosmological constant, £Ia = 0.75, the 
normalization of density fluctuations, as = 0.9 and Hubble 
constant, h = H /(l OOkms^Mpc- 1 ) = 0.73. 

Due to memory restrictions, the Fourier mesh used to 
set up the initial particle displacements has a dimension of 
1580 3 grid points which is not commensurate with the cube 
root of the particle number mesh. We therefore avoided us- 
ing a regular particle grid to set up the initial conditions, as 
this would have led to a spurious feature in the power spec- 
trum of the initial conditions at the beat frequency between 
the particle grid and th e Fourier me s h. Instead, we use d 
a glass-like distribution (jWhitei [l99~i iBaugh et all Il995l ). 
The input power spectrum of density fluctuations in linear 
perturbation theor y is calculated using the CAMB package of 
I Lewis et al.l (|2000l ). The amplitude of the Fourier modes is 
drawn from a Rayliegh distribution with mean equal to the 
linear theory power spectrum and the phase is drawn at 
random from the interval to 2-7T. The initial density field is 
generated by perturbing particles from the glass-like distri- 
bution, using the approximation of Zel'dovich (1970). 

The simulations were started at a redshift of z = 
63. The Zel'dovich (1970) approximation used to set up 
the initial pattern of density fluctuations produces tran- 
sients which can be seen in clustering signal measured 
for the dark m atter at expansion factors close to the 
starting redshift (jEfstathiou et aL 1 1 19851 : IBaugh et al.lll995l : 
ICrocce et al.l [20061 ). Later on, we will use the power spec- 
trum from a high redshift output from the simulation, 
z = 15, as a proxy for linear perturbation theory, so it is 
important to check that this power spectrum in particular, 
and also the power spectra measured at all subsequent out- 
puts are insensitive to the choice of starting redshift. We 
test this by comparing the power spectrum of the dark mat- 
ter at z = 15 in our standard run with the spectrum mea- 
sured in a test run which started at z = 127, but which 
did not run all the way through to z = 0. The top panel 
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of Fig. \T\ shows that the power spectra measured for the 
dark matter in these two cases, divided by the power spec- 
trum predicted by linear perturbation theory at z = 15. 
The fluctuations in the measured power at low wavenum- 
bers around the linear theory prediction reflect the sample 
variance noise which is not negligible even in a simulation of 
the volume of the BASICC. The lower panel in Fig. [T] shows 
the z — 15 power spectrum measured from the run started 
at 2 = 63 divided by that measured from the run started at 
2 = 127. At large wavenumbers, the effect of transients is 
visible, although quite small, ~ 1%. The focus of this paper, 
however, is the form of the power spectrum over wavenum- 
bers smaller than k = 0.4 foMpc -1 , for which the spectra 
measured at z — 15 for the two different choices of starting 
redshift agree to better than 0.3%. Our results are therefore 
unaffected by any transients resulting from the use of the 
Zel'dovich approximation. 



2.1.1 The high resolution simulation: the BASICC 

The BASICC simulation covers a comoving cubical region 
of side 1340 h~ x Mpc, in which the dark matter is rep- 
resented by more than 3 billion (1448 3 ) particles. The 
equivalent Plummer softening length in the gravitational 
force is e = 50/i _1 kpc, giving a dynamic range in length 
of almost 27,000. The volume of the computational box, 
2.41 h~ 3 Gpc 3 , is alm ost twenty times the v olume of the Mil- 
lennium Simulation (|Springel et al.l [20051 ). and more than 
three times the volume of the catalogue of luminous red 
galaxies from th e SPS S used to detect the acoustic peak by 
lEisenstein et alj (120051 ). The BASICC volume is within a fac- 
tor of two of that proposed for a survey with WFMOS at 
2 ~ 1 (jGlazebrook fc BlakdliooEh . The simulation occupied 
the full 0.5 Terabytes of RAM of the second upgrade of the 
Cosmology Machine at Durham. The run took 11 CPU days 
on 506 processors, the equivalent of 130,000 cpu-hours. 

The particle mass in the BASICC simulation is m p = 
5.49 x 10 10 /i -1 Mq. This is approximately 64 times larger 
than the particle mass used in the Millennium Simulation. 
The mass resolution limits the usefulness of dark matter 
halo merger trees from the BASICC, so we have chosen to 
output at a modest selection of redshifts: z~0, 0.3, 0.5, 1, 
2, 3, 4, 6, 8, 10, 15 and 63. Each of these outputs occu- 
pies ~ 100 Gb of disk space. : In each snapshot we have 
identified groups of dark matter par ticles using a friends-of- 
friends algorithm (|Davis et al.ll 19851 ) with a linking length of 
0.2 times the mean inter-particle separation. We have stored 
groups with 10 or more particles, i.e. haloes more massive 
than 5.49 x 10 11 h' 1 Mq. There are 17 258 579 haloes in the 
2 = output of the simulation with ten or more particles. 
The most massive halo has a mass of 6.74 x 10 15 ft -1 Mq 
and 860 haloes have a mass in excess of the Coma cluster 
(« 10 15 ft- 1 M Q ). 

The BASICC simulation sits betw een the Millennium and 
Hubble Volume (|Evrard et alj|2002h simulations. Its unique 
combination of mass resolution and volume makes it ideal for 
studying the large scale distribution of galaxies and clusters 
alike. 



2.1.2 The ensemble of low resolution simulations: 
L-BASICC 

We also generated an ensemble of 50 "low-resolution" sim- 
ulations to study the sample variance in the BASICC and to 
test an analytic model for the errors expected on measure- 
ments of the power spectrum, which we discuss in the next 
subsection. These low resolution runs (L-BASICC) have ex- 
actly the same cosmological parameters as the BASICC and 
the same box size (see Table 1), but they have fewer par- 
ticles (448 3 ). For each realization, a different random seed 
is used to set up the initial density field. The starting red- 
shift of these simulations is z — 63. The particle mass is 
compa rable to that emplo yed in the Hubble Volume simu- 
lation l|Evrard et alj|2002h . Each L-BASICC simulation took 
0.8 days to run on 16 processors of the third upgrade of 
the Cosmology Machine. The total volume of the ensemble 
is 120 h~ A Gpc 3 , more than four times that of the Hubble 
Volume, making this a unique resource for studying the fre- 
quency of rare objects in a ACDM universe. For L-BASICC, 
the position and velocity are stored for every particle at 4 
output times (z = 0.0, 0.5, 0.9, 3.8); we also produce a halo 
catalogue at each redshift retaining objects with ten or more 
particles (corresponding to a mass of 1.8 xlO 13 ft _1 M ). As 
we shall see in later sections, the ensemble allows us to assess 
whether or not a particular result is robust or simply due to 
sampling fluctuations. Due to their limited mass resolution, 
it is not feasible to populate these simulations with galaxies 
using the method outlined below (§2.2). 

2.1.3 Power spectrum estimation and errors 

The two point statistics of clustering, the correlation func- 
tion, and its Fourier transform, the power spectrum, P(k), 
are the most commonly employed measurements of clus- 
tering. In this paper we focus on the power spectrum; in 
Sanchez et al. (2007, in preparation), we address the vis- 
ibility of the acoustic oscillations in the correlation func- 
tion. The standard way to quantify the amplitude of a 
density fluctuation is by means of the density contrast, 
6(x,t) = (p(x,t) — p)/p. If we consider the Fourier trans- 
form of the density contrast, p^, then the power spectrum 
is defined as the modulus squared of the mode amplitude, 

P(k) = (|5 k )| 2 ). 

There are two steps in the computation of the power 
spectrum from a distribution of discrete objects, such as 
dark matter particles, dark haloes or galaxies. Firstly, a den- 
sity field is constructed by assigning the objects to mesh 
points on a cubic grid. In the simplest mass assignment 
scheme, the nearest grid point, the contribution of each ob- 
ject to the density field is confined to the cell in which it 
is located. In higher-order assignment schemes, the mass of 
the particle is shared with adjacent cells. Here, we use the 
cloud -in-cell assignment scheme (see iHocknev fc Eastwood! 
Il98ll) . Secondly, we perform a Fast Fourier Transform of the 
density field. The power spectrum is obtained by spherically 
averaging the resulting Fourier mode amplitudes in annuli 
of radius 5k = 2n/L = 0.0047 h Mpc -1 . 

The mesh we use to store the density field has Np FT = 
512 3 grid points. Estimating the density on a grid al- 
ters the form of the power spectrum at wavenumbers ap- 
proaching the Nyquist frequency of the grid (fcNyquist = 
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Figure 2. The power spectrum of the dark matter in real-space 
measured at the starting redshift of the BASICC, z = 63 (red 
points). The corresponding prediction of linear perturbation the- 
ory is shown by the green (solid) line. The blue (dot-dashed) curve 
shows the power spectrum of the unperturbed glass-like distribu- 
tion of particle positions. The dashed line shows the Poisson noise 
expected for the number density of dark matter particles used in 
the BASICC. The noise of the initial particle distribution is much 
less than Poisson. The arrow marks the position of the Nyquist 
frequency of the FFT grid. 



2n/L TVfft/2 = 1.2/iMpc" 1 in our case). The degree of 
modification and the precise wavenumber above which the 
power spectrum is distorted de pend upon the choice of as- 
signment scheme (Hatton 1999: Ijingll2005h . In practice, for 
the size of FFT mesh we use, this has little impact on the 
recovered power spectrum for wavenumbers of interest; the 
measured amplitude differs by less than 1% from the true 
value at a wavenumber k ~ 0.8/iMpc -1 ; in most cases we 
focus on the form of the power spectrum on large scales, 
k < 0.4 h Mpc -1 . Nevertheless, we correct for the effects of 
the cloud-in-cell assignment scheme by dividing each mode 
by the Fourier transform of a cubical top hat: 



3(k x , ky , kz^j 
where 



&\kx j ky , kz*) 



sinc( 



fc x L 
2iVpp"p 



lsinc (2^|7) smc ( 



k z L - 
2iV FFT - 



(1) 



sinc(:r) = 



sin(:r) 



(2) 

Note this is different from the approach taken by Ijind 
(2005), who applied a correction to the spherically averaged 
power spectrum. 

A further possible distortion to the form of the mea- 
sured power spectrum is discreteness noise and the asso- 
ciated Poisson or shot noise. Poisson-sampling a continu- 
ous density field with point objects of space density, n, in- 
troduces a spurious contribution that should be subtracted 
from the measured power spectrum: P CO rr(fc) = Pmeas(fc) — 
1/n. In the case of dark matter halo centres or galaxies, the 
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Figure 3. The fractional error in the power spectrum of the 
dark matter (top panel) and in the power spectrum of haloes 
more massive than 1.8 X 10 13 h~ 1 Mq (bottom panel), estimated 
using the low resolution simulations from the dispersion of P(k) 
around the ensemble mean. The smooth black curves show the er- 
ror predicted by the analytical expression given in Eq. \3\ The red 
points show the scatter from the ensemble of low resolution sim- 
ulations. The arrow in the bottom panel shows the wavenumber 
for which nP(k = O^hMpc -1 ) = 1. 



need for such a correction is justified. However, in the case of 
dark matter particles in our simulations, one should not sub- 
tract Poisson shot noise from the power spectrum because 
the particles were initially laid down by perturbing a glass- 
like configuration which is sub-Poissonian in nature. This is 
clear from Fig. [2] which shows the power spectrum measured 
for the dark matter in the initial conditions of the BASICC. 
The red curve shows the spectrum measured in the simula- 
tion and the smooth green curve shows the input spectrum 
predicted by linear perturbation theory. The two agree re- 
markably well over a wide range of wavenumbers. The power 
spectrum of the unperturbed glass-like particle distribution 
is shown by the blue curve. For the wavenumbers of interest, 
the power spectrum of the glass is many orders of magnitude 
below the discreteness noise expected for a Poisson distri- 
bution of objects with the same space density as the dark 
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matter particles, as shown by the dashed line. In this paper, 
we do not apply any shot noise correction to power spec- 
tra measured for the dark matter, but we do make such a 
correction for spectra estimated for samples of haloes and 
galaxies. 

To close this subsection, we turn our attention to the 
error on the measurement of the power spectrum. A com- 
monly used expression for the fractional error in the mea- 
sured power spectrum was derived by Feldman, Kaiser & 
Peacok(1994) (see also Efstathiou 1988, for a similar argu- 
ment applied to the two point correlation function): 



a 
P 



^modcs 



(3) 



where n mo d cs is the number of Fourier modes present in a 
spherical shell of width 8k, which depends upon the sur- 
vey volume V: for k » 2tt/V 1/3 , this is given by 

^modcs 

V4nk 2 5k/ (2n) 3 . The first term on the right hand side of 
Eq. [3] quantifies the sample variance in the measurement, 
which decreases as the square root of the number of modes 
or, equivalently, as the square root of the volume probed. 
The second term arises from the discreteness of the ob- 
jects under consideration. The combination Pn quantifies 
the amplitude of the power spectrum in units of the Poisson 
shot noise, effectively giving the contrast of the power spec- 
trum signal relative to the shot noise level. In the case where 
Pn S> 1, ct/P oc 1/k. On the other hand, when the ampli- 
tude of the power spectrum is comparable to the shot noise, 
and if P(k) oc A; -1 , then the fractional error in the power is 
approximately independent of wavenumber. We have tested 
this prescription in both regimes against the diagonal ele- 
ment of the covariance between power spectrum measure- 
ments extracted from the ensemble of low resolution sim- 
ulations, as shown in Fig. [3] Over the wavenumber range 
of interest, the agreement is reasonably good for samples in 
which the shot noise is negligible compared to the cluster- 
ing signal. For samples with low contrast power measure- 
ments, such as is the case for dark matter haloes used in 
the bottom panel of Fig. [3] the analytic expression works 
well until k ~ 0.1/iMpc -1 and then overpredicts the errors 
by up to 50%. We note that nonlinearities and the impact 
of the window function of a realistic survey could introduce 
off-diagonal terms in the power spectrum covariance matrix. 
In Section 5.3, we compare the constraints on the recovered 
oscillation scale using the scatter from the ensemble and 
using the simple mode-counting argument outlined above. 
We find good agreement which suggests that mode-coupling 
does not make a significant contributions to the errors on 
the scales relevant to the BAO. 
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Figure 4. Upper panel: the fraction of 'resolved galaxies' in the 
high resolution N-body simulation as a function of magnitude, at 
different output redshifts (as given by the key in the lower panel) . 
The magnitude is in the observer-frame H-band; to obtain an ap- 
parent _R-band magnitude, the distance modulus corresponding to 
the redshift should be added to the plotted magnitude. The verti- 
cal lines mark the magnitude at which the galaxy sample is 100% 
complete at each redshift. Lower panel: the cumulative luminos- 
ity function of galaxies brighter than a given i?-band magnitude, 
for different redshifts as given in the key. The vertical lines show 
the 100% completeness limits at each redshift and the horizontal 
lines indicate the associated space density of galaxies. 



2.2 Modelling the formation and evolution of 
galaxies 

The N-body simulations described in the previous section 
follow the growth of fluctuations in the mass which is dom- 
inated by collisionless matter. To connect the predictions of 
the cold dark matter theory to forthcoming galaxy surveys, 
we need to predict which structures host galaxies and how 
galaxy properties depend on halo mass. 

Some authors have chosen to incorporate galaxies into 
an N-body simulation empirically by using a parametric 



model called a halo occupation distribution function (HOD) 
to describe the probability distribution of galax ies expected 
in haloes of a given mass (|Benson et al . 2000). The form 
of the HOD is constrained to reproduce a particular clus- 
teri ng measurement, such as the galaxy correlation function 
(e.g. jPeacock fc Smithl l2000l; ISeliakll200ol ; IScoccimarro et all 
l200ll ; ICoorav fe Shethl I2002T ) . This approach has been ap- 
plied to the study o f the detectability of acoustic oscillations 
by several authors (ISeo fc Eisensteinll2005l : ISchulz fc White! 
2006;|Huff et al ] |2006l 1. Two assumptions are made when us- 
ing the HOD to populate an N-body simulation with galax- 
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ies. Firstly, the parameterization used for the HOD is as- 
sumed to provide an accurate description of the manner in 
which galaxies populate haloes across a wide range of halo 
mass. Detailed comparisons between the clustering predic- 
tions made using HODs and those obtained directly from 
simulations of galaxy form ation show that in practice, the 
HOD s do a reasonable job (jBerlind et a l, 2003; Zhcn get al.l 
2005). Recently, one of the fundamental assumptions which 
underpins the HOD approach has been called into q uestion. 
Using the Millennium simulations. iGao et al.l (|2005h demon- 
strated that the clustering of dark matter haloes depends 
on a second parameter, such as the for mation time of the 
halo, in addition to halo mass (see also Harker et al.ll2006l 



and lWechsler et~al ]2006t I Wetzel et al.ll2007l ). In practice, for 
typical galaxy samples, this effect is lar gely washed out du e 
to the mix of halo properties sampled ( Croton et al .1120071 ). 
The second implicit assumption in the HOD method when 
applied to an N-body simulation is that all of the haloes 
in which galaxies are expected to be found can be resolved 
in the simulation; if the mass resolution of the simulation 
turns out to be inadequate, then the HOD realized will be 
distorted to compensate, compared with the true, underly- 
ing HOD in the Universe. 

In this paper, we take a more physical approach and 
make an ab initio prediction of which dark matter haloes 
should contain galaxies by modelling the physics of the bary- 
onic component of the universe. We do this using a semi- 
analytic model of galax y formation (for a review of this tech- 
nique see lBaughll2006h . The semi-analytic model describes 
the key physical processes which are thought to determine 
the formation and evolution of galaxi es. We use the GALFORM 
code introduce d by Cole et al. 1120001) and developed in a se- 
ries of papers dBenson et al.ll2002l . 120031 ; iBaugh et al1l2005l ; 
I Bower et al. 20061). The speci fic model we use is the one pro- 
posed by Baugh et a l. (2005), which reproduces the abun- 
dance of Lyman-break galaxies at z = 3 and z = 4, the 
number counts of sub-mm detected galaxies (with a median 
redshift z ~ 2), and a rough match to the abundance of lu- 
minous red galaxies (Almeida et al. 2007, in preparation), 
whilst at the same time giving a reasona ble match to the 
observ ed properties of local galaxies (e.g. iNagashima et all 
l2005al lbl; I Almeida et ai1l2007l ). 

A key advantage of using a semi-analytic model is that 
we can investigate how the manner in which galaxies are 
selected affects the accuracy with which the acoustic oscilla- 
tions can be measured. The model predicts the star forma- 
tion history of each galaxy and uses this to compute a spec- 
trum, broadband magnitudes and emission line strength s 
(for examples of the latter, see lLe Delliou et alj |2005. 2006). 
We can therefore select samples of model galaxies by apply- 
ing precisely the same criteria which will be applied in the 
proposed surveys. 

Our methodology mirrors t he h ybrid schemes intro - 
duced bv lKauffmann et all i|l997l ) and lBenson et al.l l|2000l) . 
We use a Monte Carlo technique to generate merger trees for 
dark mater haloes since our simulation outputs do not have 
the resolution in time or mass necessary to allow the con- 
struction of merger trees. (Sec Baugh 2006 for a discussion 
of the relative merits of these two approaches.) 

We first construct a grid of halo masses at the redshift 
of interest, which extends to lower mass haloes than can 
be resolved in the simulation. We then generate a number 



of Monte-Carlo realizations of mass assembly histories for 
each mass on th e grid, using the algorithm introduced by 
ICole et al l l|2000l) . The number of realizations is chosen to 
allow robust predictions to be made for observables such as 
the galaxy luminosity function. The halo merger history is 
input into the semi-analytic code and the properties of the 
galaxy population are output at the redshift for which the 
galaxy catalogue is to be constructed. In the calculations in 
this paper, we output the broadband magnitudes in the R, I 
and K bands and the equivalent widths of H a and Oil [3727] 
for each galaxy. Finally, haloes from the grid are matched 
with haloes of similar mass identified in the N-body simula- 
tion. The central galaxy in each halo is assigned to the centre 
of mass of the matched halo in the simulation. The satellite 
galaxies are assigned randomly to dark matter particles in 
the halo. Galaxies placed in the simulation box in this way 
are called 'resolved galaxies'. The Monte Carlo merger trees 
will not, of course, correspond in detail with those of the 
matched halos in the N-body simulation. Howe ver, to the 
extent that the halo assembly bias discussed bv lGao et alj 
(2005) can be neglected, the properties of the trees are sta- 
tistically similar for haloes in the same mass range. 

Because of the finite mass resolution of the N-body sim- 
ulation, galaxy samples generated by populating resolved 
haloes will be incomplete fainter than some magnitude limit. 
In principle, since we are using Monte-Carlo merger trees, 
we can follow galaxies down to arbitrarily faint magnitudes 
within a resolved dark matter halo. However, as we con- 
sider progressively fainter objects, some fraction of these 
galaxies should also appear in haloes which the simulation 
cannot resolve, causing the sample to become incomplete. 
Thus, in some instances we need to consider galaxies which 
we would expect to find in haloes below the mass resolu- 
tion of the simulation. These galaxies are called "unresolved 
galaxies" and are placed in the box in the following way. 
A volume-limited sample of galaxies is generated using the 
semi-analytic model, with a volume equal to that of the sim- 
ulation cube. Only galaxies which reside in haloes from the 
grid which are less massive than the resolution limit of the N- 
body simulation are considered. (Recall that the grid of halo 
masses used in the semi-analytic calculation extends to lower 
mass than those resolved in the simulation). These galax- 
ies are assigned to randomly selected dark matter particles 
which have not been identified as members of halos identi- 
fied by the friends-of-friends algorithm. This ap proach was 
adopt ed for one of the mock catalogues used in I Cole et all 
( 200|). As we will see below, the unresolved galaxies are a 
minority within any of the samples we consider. They have 
little effect on the measured power spectrum, producing only 
a modest change in the amplitude of the clustering signal. 

We can use the semi-analytic calculation carried out on 
the grid of halo masses to find the completeness limit of the 
galaxy catalogue in the N-body simulation. To do this, we 
use the galaxy formation calculation carried out using the 
grid of halo masses to compute the cumulative luminosity 
function of galaxies, starting with the brightest galaxy, for 
two cases: 1) without any restriction on the mass of the 
halo which hosts the galaxy and 2) considering only those 
galaxies which reside in haloes above the resolution limit 
of the simulation. We then divide the second estimate of 
the cumulative luminosity function by the estimate made 
without any restriction on halo mass. 



8 Angulo et al. 



10 5 




irr 2 l i i i i 

0.01 0.10 1.00 

k [h Mpc -1 ] 



Figure 5. The growth of the power spectrum of density fluctua- 
tions in the dark matter, as measured in real-space. The smooth 
curves show the predictions of linear perturbation theory at the 
redshifts indicated by the key. The power spectra measured in the 
low resolution ensemble at z = are plotted to show the sampling 
variance for a simulation box of side 1340 /i -1 Mpc. The smallest 
wavenumber plotted corresponds to the fundamental mode in the 
simulation, 2tt/L = 0.0469 h' 1 Mpc. The maximum wavenumber 
shown is 0.67 times the Nyquist frequency of the FFT grid, chosen 
to avoid any aliasing effects. 



The completeness ratios calculated in this way are 
shown for z = 0, 1 and 2 in Fig. [4] The vertical lines show 
the magnitude limit down to which the 'resolved galaxy' 
catalogues are 100% complete. The lower panel shows the 
cumulative luminosity function in the model at the same 
redshifts, with horizontal lines marking the space density of 
galaxies at the sample completeness limit. (The magnitudes 
plotted are observer-frame absolute magnitudes in the R- 
band. The apparent magnitude is obtained by adding the 
appropriate distance modulus for each redshift. All magni- 
tudes are on the AB scale.) The 2 = 2 sample is complete 
down to Mr — 5 log h — —23, or, equivalently to a space 
density of 3.2 x 10 -5 h 3 Mpc -3 . Faintwards of this magni- 
tude, the completeness drops sharply to around 30 — 40%. 
The situation is much more encouraging at z — 1. Here, the 
galaxy catalogue is complete to Mr — 5 log h = —22.3 (cor- 
responding to a space density of just under 10 -4 ft 3 Mpc -3 ) 
and faintwards of this there is a much more modest drop 
in the fraction of galaxies resolved in the simulation. The 
simulation resolves around two thirds of the space density 
of galaxies expected in the proposed WFMOS survey. At 
2 = 0, the galaxy samples are complete to a much higher 
space density, in excess of 10 -3 /i 3 Mpc -3 . 



3 THE POWER SPECTRUM OF GALAXY 
CLUSTERING 

In this section we examine the various phenomena which are 
responsible for changing the form of the power spectrum of 
galaxy clustering from that expected in linear perturbation 
theory. We systematically add in new effects and elements of 
sample selection, considering first the power spectrum of the 
dark matter, looking at nonlinear evolution ( i|3.1|l and the 
impact of peculiar velocities f £!3.2[) . before moving onto dark 
matter haloes f £!3.3[) and finally to synthetic galaxy samples 

(MM- 

For completeness, we first explain some of the terminol- 
ogy we use in this section. There are three types of phenom- 
ena responsible for distorting the linear theory power spec- 
trum: i) non-linear growth of fluctuations, ii) redshift-space 
distortions and iii) bias. Non-linear growth refers to the cou- 
pled evolution of density fluctuations on different scales. 
Redshift-space distortions describe the impact of gravita- 
tionally induced peculiar motions on the clustering pattern. 
We will refer to clustering measurements as being made in 
"real-space" or "redshift-space"; in the latter case peculiar 
motions are taken into account, as we describe in §3.2. The 
term "bias" has a range of meanings in the literature. Bias 
is used to describe the boost in the clustering of a particular 
tracer (e.g. galaxies or clusters) relative to a reference point, 
which could be the clustering of the dark matter in either 
linear perturbation theory or taking into account nonlinear 
evolution. One of the earliest uses of the concept of bias was 
in the application of the high peaks model to exp lain the 
enhanced clustering of Abell clusters (Kaiser 198J). In this 
model, clusters are associated with rare peaks in the initial, 
Gaussian density field. The bias is defined as the square root 
of the ratio of the two-point correlation function of peaks 
of a certain minimum height to the clustering of the mass 
expected in linear perturbation theory. When considering 
galaxies, it is perhaps more natural to think in terms of a 
modulation of clustering relative to that displayed by the 
underlying mass at the same epoch, since galaxies populate 
dark matter haloes. In this case, the galaxy clustering will be 
measured relative to that of the evolved matter distribution. 
On large scales, these two reference points, the clustering of 
the matter expected in linear perturbation theory or the 
evolved clustering, should be essentially the same. We shall 
see later that this is approximately the case for the scales 
over which we compare clustering signals to measure bias 
factors. 



3.1 The nonlinear growth of matter fluctuations 

The early stages of the growth of a density fluctuation are 
particularly simple to describe analytically. The fluid equa- 
tions can be written in terms of the perturbation to the 
density and Fourier transformed. In the simplest case, when 
the density contrast 5 <C 1, the Fourier modes evolve inde- 
pendently of one another. This is called linear growth. In 
this regime, the power spectrum changes in amplitude with 
time, but not in shape. The shift in amplitude is described 
by the growth factor D, which is a function of the densities 
of matter and dark energy (as quantified by the present day 
density parameters, Qm and Q a, for matter and dark ener gy 
respectively) and redshift (see Ifleathll 19771 ; |Peebleslll980t ): 
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Figure 6. The nonlinear growth of the power spectrum. Here we 
divide the power spectrum in real-space measured at the redshift 
indicated by the key by the power spectrum at Z. = 15, after taking 
into account the change in the growth factor. Any deviation of 
the resulting ratio from unity indicates a departure from linear 
perturbation theory. The dashed lines show the same ratio as 
predicted using the ansatz of Smith et al. (2003). 



P(k, z) = D 2 (z, n M ,^A)P(k, z = Q), 



(4) 



where D(z = 0) = 1. 

We plot the power spectrum of the dark matter in real- 
space measured from the BASICC at different output redshifts 
in Fig. [S] The approximately linear growth of the power 
spectrum is readily apparent on large scales (low k). In an 
Einstein - de Sitter universe (Qm = 1), the growth factor is 
equal to the expansion factor. If dark energy plays a role in 
setting the rate at which the universe expands, the growth 
of fluctuations is suppressed relative to the Einstein - de Sit- 
ter case at late times. The BASICC started at z B = 63, so if 
JIm = 1, we would expect to see the power spectrum grow 
in amplitude by a factor of (1 + z B ) 2 = 4096 by z = . Usin g 
the approximate formula provided bv lCarroll et al. I (|l992h . 
we expect a suppression in the growth of the power by a 
factor of 0.5537 for the cosmological parameters used in the 
simulation. This gives an overall growth in power from the 
initial conditions to the present of a factor of 2268. This 
agrees to within 0.6% with the factor expected from a di- 
rect numerical integration of the equation givin g the growth 
factor (eqns. 28 and 9 from lCarroll et alii 19921 ). which gives 
2281.01. In the simulation, we find that the power in the 
fundamental mode grows by a factor of 2285.21 from the 
initial conditions at z = 63 to z = 0, which agrees with the 
growth predicted by linear perturbation theory to 0.02%. 

Fig. [S] shows that the growth of the power spectrum 
is clearly not linear at high wavenumbers. The shape of the 
spectrum at high k at late times is different from that at high 
redshift, because the growth of modes of different k becomes 
coupled. This behaviour can be followed to some ex tent us- 
ing second- and higher-order perturbation theory |Peeblesl 



trast approach es unity, second-order p erturbation theory 
breaks down (IBaugh fc Efstathioulll994r ). The coupled evo- 
lution of the Fourier modes starts on surprisingly large 
scales, which demonstrates the necessity of a large volume 
simulation to accurately follow the development of the power 
spectrum (Smith, Scoccimarro & Sheth 2006). This can be 
seen more clearly if we divide the measured spectrum by 
the growth expected according to linear perturbation the- 
ory, as is done approximately in Fig. [6] In this plot, we 
have divided the power spectra measured from the simula- 
tion by the spectrum measured at z = 15, scaled by the 
square of the appropriate growth factor. This reduces the 
noise in the ratio arising from the finite number of modes 
realized at small wavenumbers in the simulati on volume 
jBaugh fc Efstathioul 1 1994 ISpringel et"aH 120051 ). Any de- 
viation away from unity signifies a departure from linear 
perturbation theory due to coupling between modes. The 
ratio shows a characteristic dip at low k, i.e. less power than 
expected in linear theory, bef ore showing a strong enhan ce- 
ment at higher wavenumbers l|Baugh fc Efstathioull 199*4 ) . It 
is remarkable that the transition between a deficit and excess 
of power happens at the same wavenumber, k ~ O.lfoMpc -1 , 
at different epochs. The suppression in power at low k, on the 
order of a 3%, is not as strong as that seen in an Einstein - 
de Sitter universe (see figure 4 of lBaugh fc Efstathiou|[l994 ) . 
Nevertheless, this drives the spectacular boost in power seen 
at higher wavenumbers. The dip in power is largest around 
k ~ 0.05ftMpc -1 , which corresponds to a length scale of 
2%/k ~ 125 h" 1 Mpc, close to the wavelength of the acoustic 
oscillations. Several authors have proposed ansatzes which 
transform the linear perturbation theory power spectrum 
into the non- linear power spectrum (e.g. lHamilton et al.l 
1 199 it iPeacock fc Dodddll994 Il996l : ISmith et alj|2003h. We 
plot t he predictions of the model proposed bv lSmith et al.l 
( 2003) in Fig. © using dashed lines. The ratio is computed 
by dividing the power spectrum at the epoch of interest by 
the suitably scaled prediction of the model for z=15. The 
agreement is excel lent at high r edshif t. At z — 0, at higher 
wavenumbers, the ISmith et al.l l|2003l ) formula recovers the 
simulation results to within 5% over the range plotted. 



on 



3.2 The impact of redshift-space distortions 
the power spectrum 

In a spectroscopic galaxy survey, the radial distance to an 
object is inferred from its measured redshift. The shift in 
the spectral features of the galaxy is produced by two con- 
tributions to its the apparent velocity: the expansion of the 
universe, which is responsible for the Hubble flow at the 
true distance to the galaxy, and local inhomogeneities in the 
gravitational field around the object, which generate an ad- 
ditional, "peculiar" velocity. Since we cannot correct a priori 
for the effects of the local gravitational field when inferring 
the radial distance from the Hubble law and the measured 
redshift, an error is made in the distance determination. The 
impact of such errors on the form of the measured power 
spectrum of clustering is called the redshift-space distortion. 

Peculiar motions display two extremes which produce 
different types of distortion to the power spectrum: i) On 
large scales, coherent bulk flows out of voids and into over- 
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Figure 7. The ratio of the power spectrum measured for the 
dark matter in rcdshift-space, i.e. including the impact of pecu- 
liar motions in the distance determination, to the power spec- 
trum measured in real-space. The deviation from unity shows the 
redshift-space distortion to the nonlinear power spectrum. The 
results are shown for selected output redshifts, as indicated by 
the key. The horizontal dotted lines indicate the boost in the 
redshift-space power expected due to coherent flows, as predicted 
by Eq. \E\ The dashed lines show a simple fit to the distortions 
(see eq(6)l. 



dense regions lead to an enhancement in the density inferred 
in red shift-sp a ce, a nd hence to a boost in the recovered 
power. iKaiserl (|l987l ) derived a formula for the enhancement 
of the spherically averaged power, under the assumption of 
linear perturbation theory for an observer situated at infin- 
ity (the plane parallel approximation): 



/ 



Ps(k) 
Pr(k) 



(l + |/3+i/3 2 ) 



(•») 



where P s (k) is the power spectrum in redshift-space, P T (k) 
is the spectrum in real-space and (5 — (dlog 8/d log a) /b ~ 
Q l ^(z)/b, where b is the bias factor (6=1 for the dark mat- 
ter; for a discussion of the dependence of the growth factor 
on JIm, see Linder 2005; Linder & Cahn 2007). ii) On small 
scales, the random motions of objects inside virialized dark 
matter haloes cause structures to appear elongated when 
viewed in redshift-space. leading to a damping of the power. 
iPeacock fc Doddsl l|l994l ) discussed a model for the redshift- 
space power spectrum, whic h takes into accoun t both limits 
of peculiar motions (see also IScoccimarr ol l2004h . 

Fig. [7] shows the ratio of the power spectrum measured 
for the dark matter in redshift-space to that measured in 
real-space, at redshifts z = 3, 1 and 0. The dotted lines in- 
dicate the boost expected in the red shift-space p ower, com- 
puted using the expression in Eq. [5l (|Kaiserlll9"87T ). This fac- 
tor changes with redshift because the matter density param- 
eter is changing. Fig.[7]shows that this behaviour is only ap- 
proached asymptotically, on scales in excess of 100 ft" 1 Mpc. 



At higher wavenumbers, the power measured in redshift- 
space is suppressed by random motions. The dashed lines in 
this plot show a simple fit to this ratio 



/ = 



(6) 



where a is a free parameter, which is loosely connected to the 
pairwise velocity dispersion. The degree of damping grows 
between z = 3 and z = 1, but changes relatively little by 
z = 0. We shall see in later sections that the form of the 
redshift-space distortion to the power spectrum depends on 
the type of object under consideration. 



3.3 The power spectrum of dark matter halos in 
real and redshift-space 

In modern theories of galaxy formation, dark matter haloes 
play host to galaxies. It is therefore instructive to compare 
the power spectra measured for different samples of haloes 
to that of the dark matter as a step towards understanding 
the power spectrum of galaxies. 

A common conception is that the clustering of haloes 
is a scaled version of the clustering of the underlying mass, 
with the shift in clustering amplitude quan tified in terms 
of a bias factor, b, wher e b 2 = Phaios/^dm l|Cole fc Kaiserl 
1 19891 ; iMo fc White! Il996h . As we commented earlier, since 
we use the dark matter power spectrum on large scales to 
define a bias, this is approximately the same as using the 
linear perturbation theory spectrum. Many authors have 
tested analytical prescriptions for comp uting the bias pa- 
ramet er using extensions of the theory of lPress fc Schechterl 
dl974) (e.g. Mo Jing fc White 1997; Sheth, Mo fc Tormen 
2000; |jindil99Sl; iGovernato et al.|[l999l ; IColberg et aT1l2000l ; 
ISeliak fc Warrenl I2004T ). In the extended Press-Schechter 
theory, the bias is only a function of halo mass and redshift. 
However, recent analyses of high resolution, large volume 
simulations have revealed some dependence of halo cluster- 
ing on a second parameter besides mass, such as the halo's 
formation redshift or concentration paramet er (|Gao et alj 
120051 ; lHarker et aT1l2006l ; IWechsler et~al]|2006l '). 

In Fig. [8] we show that this simple picture, in which the 
clustering of haloes is a shifted version of that of the dark 
matter, is actually a poor approximation to what we find 
in the simulation. We show the ratio of the power spectrum 
of a sample of dark matter haloes measured in real-space 
to a scaled version of the linear perturbation theory power 
spectrum. The amplitude of the linear theory spectrum used 
in the ratio takes into account the growth factor appropri- 
ate to the output redshift and an effective bias, which is set 
by matching the linear theory prediction for the mass spec- 
trum to the measured halo spectrum on large scales, i.e. 
for wavenumbers in the range 0.0046 < (fc/hMpc- 1 ) < 0.1. 
Each panel in Fig. [8] corresponds to a different output red- 
shift from the simulation. For each redshift, we have defined 
three samples of dark matter haloes, which contain the same 
number of objects. The mass intervals are set relative to the 
average halo mass present in the respective outputs, with 
"low", "mean" and "high" mass samples considered. Each 
of these contains 20% of the total number of haloes present 
at each epoch, with the mass ranges used at each redshift 
indicated on the keys. The effective bias factors of the halo 
samples are also written in the key. For comparison, the 
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Figure 8. The power spectrum of dark matter haloes measured in real-space compared to a scaled version of the prediction of linear 
perturbation theory, which takes into account the growth factor and an effective bias computed on large scales k < 0.1/iMpc - . Each 
panel corresponds to a different output redshift. Different mass samples are considered, as indicated by the key, which correspond to 
low, average and high masses, defined in terms of the average halo mass present at each output time. The black dashed line shows the 
real-space power spectrum of the mass divided by the appropriate linear perturbation theory prediction. 




Figure 9. The power spectrum of dark matter haloes measured in redshift-space divided by the power spectrum measured in real-space 
for the same sample. Each panel corresponds to a different output redshift. Different mass samples are considered, as indicated by the 
key, which correspond to low, average and high masses, defined in terms of the average halo mass present at each output time. The 
horizontal dotted lines show the expected ratio for the boost in the amplitude of the redshift-space power spectrum due to coherent 
flows, computed using an effective bias factor estimated on large scales. The dashed lines show the best fit model of Eq. [6] which turns 
out to be a poor description of the redshift-space distortions. No suitable fits were obtained at z = 3. 



dashed line in each panel shows the corresponding ratio for 
the dark matter. 

Fig. [5] shows that at z = 3, all of the haloes consid- 
ered have effective biases much greater than unity, indicat- 
ing they are more strongly clustered than the mass. This 
situation is reversed at z = 0. At this epoch, the halo mass 
resolution of the BASICC is smaller than the corresponding 
value of MJ3 (= 5.78 x 10 12 h' 1 M at z=0). The z = 

1 M, is a characteristic mass scale defined as the mass within a 



samples have a bias of unity or smaller. In addition to the 
difference in the effective bias parameters, the shape of the 
spectrum of the haloes in these extremes is also different 
(see also Smith et al. 2006). The plot shows the shape of the 
power spectrum, after accounting for the effective bias on 
large scales. Any difference between the curves plotted for 

sphere for which the rms variance in linear perturbation theory 
is (t(M) = <5 cr it(.z), where (5 cr ; t is the extrapolated critical linear 
overdensity given by the spherical collapse model at redshift z. 
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the haloes and that for the dark matter (dashed line) shows a 
difference in the clustering signal over and above that quan- 
tified by a constant effective bias. Similar behaviour was 
found for samples of cluster mass haloe s in the Hubble Vol- 
ume simulation by I Angulo et~aH l|2005l ). 

We now consider the clustering of haloes as viewed 
in redshift-space, taking the centre of mass velocity of the 
halo as its peculiar velocity. In Fig. |9j we plot the ratio 
of the redshift-space power spectrum for the halo samples 
used in Fig. [8] to the power spectrum measured in real- 
space. As we did before for the case of the dark matter 
(Fig. [7]), we indicate the boost in power expected on large 
scales (small k) due to coherent bulk flows of haloes. The 
boost is calculated from Eq. [5] using the effective bias of 
the halo sample. The plot shows that the redshift-space 
power spectrum at low wavenumbers is in reasonable agree- 
ment with this simple model. However, a range of behaviour 
is seen at higher wavenumbers. For haloes comparable to 
M*, the boost in power in redshift-space is less than pre- 
dicted by Eq.[5] For the more extreme, massive haloes, there 
is actually more power in redshift-space than is suggested 
by Kaise r's formula. This "exce ss" power was previousl y 
noted bv lPadilla fc Baughl (|2002r i and lAngulo et al.1 (|2005l ). 
The Kaiser formula assumes linear perturbation theory and 
breaks down in the case of objects with strongly nonlinear 
clustering. In the case of the less extreme haloes, the reduc- 
tion in power is not due to virialized motions of haloes within 
larger structures. The halo finder we have used is designed 
to return an overdensity corresponding to virialized struc- 
tures and not substructures. If the haloes were really part of 
a larger structure and were executing random motions, the 
group finder would simply have lumped them together as one 
larger structure. We are perhaps seeing instead haloes that 
have started to merge with one another, and whose motions 
have broken away from a coherent large scale flow. We know 
of no analytical description of the redshift-space clustering 
of dark matter haloes which explains this behaviour. 

3.4 The power spectrum of galaxies 

The galaxy power spectrum can be very different from the 
power spectrum of a sample of dark matter haloes. The way 
in which the galaxies are distributed among haloes changes 
the form of the power spectrum. In a mass- limited sample of 
haloes, the contribution of each halo to the power spectrum 
can be determined through its space density, which acts as 
a weighting factor when computing the contribution of the 
halo to the clustering signal. The number of galaxies per halo 
acts to modify this weight e.g. more massive haloes could 
contain more galaxies than less massive haloes. Furthermore, 
the presence of satellite galaxies within a halo means that 
one expects to see a damping in power on small scales in 
redshift-space, due to the random motions of the satellites 
within the virialized dark halo. The precise modification of 
the power spectrum depends in detail on how galaxies pop- 
ulate dark matter haloes. As we discussed in §2.2, we have 
carried out an ab initio calculation of the number of galaxies 
per halo, using a semi-analytic model of galaxy formation. 
We are able to predict observable properties of galaxies, such 
as broadband magnitudes and the strength of emission lines. 
We consider a range of galaxy samples, defined either by a 
magnitude limit alone (set in the R-band) or by combining 



an R-band magnitude limit with a colour selection (in R-I) 
or a cut on the strength of the OII[3727] emission line: 

• Sample A: magnitude-limited to reach a space density 
of 5 x I(T 4 /i 3 Mpc- 3 . 

• Sample B: magnitude- limited to reach half the space 
density of sample A, i.e. 2.5 x I0~ 4 h 3 Mpc -3 . 

• Sample C. The reddest 50% of galaxies from sample A, 
using the R — I colour. 

• Sample D. The 50% of galaxies from sample A with 
the strongest emission lines, using the equivalent width of 
Oil [3727]. 

• Sample E. The bluest 50% of galaxies from sample A, 
using the R — I colour. 

• Sample F. The 50% of galaxies from sample A with 
the weakest emission lines, using the equivalent width of 
Oil [3727]. 

The power spectra measured in real-space from the var- 
ious galaxy samples are plotted in Fig. If 01 The spectra have 
been divided by the linear perturbation theory power spec- 
trum multiplied by the square of an effective bias factor, 
which was estimated by comparing the galaxy spectra to the 
power spectrum measured for the dark matter for wavenum- 
bers k < O.f/iMpc -1 . In all cases, for the space densities we 
have chosen, the effective bias factors estimated for the sam- 
ples are modest. For comparison, the ratio of the power spec- 
trum of the dark matter in real-space to the linear theory 
prediction is also plotted, using a dashed line. The deviation 
of the dashed line from unity shows where nonlinear effects 
are important for the dark matter. Any differences between 
the plotted ratios for galaxies and mass indicate a scale de- 
pendent bias. The comparison between the dashed and solid 
curves in Fig [TU] shows that a constant bias is only a good 
approximation on large scales, k < 0.15/iMpc -1 . 

The redshift-space distortion in the galaxy power spec- 
trum is shown in Fig If f I where we plot the ratio of the 
redshift-space spectrum to the real-space spectrum for the 
galaxy samples shown in Fig. If 01 The horizontal lines show 
the Kaiser boost (Eq. expected for the effective bias of 
the galaxy sample. This ratio is only attained on the very 
largest scales and seems to be an overestimate of the size of 
the effect at z = 1. The damping of the power on intermedi- 
ate and small scales is readily apparent and, unlike the case 
with dark matter haloes, is well described by the form given 
in Eq. H 



4 CONSTRAINING THE DARK ENERGY 
EQUATION OF STATE 

In this section we outline the procedures we follow to place 
constraints on the dark energy equation of state parame- 
ter, w, by measuring the length scale imprinted by baryonic 
acoustic oscillations on the power spectrum of the various 
tracers of the density field. The transformation of a mea- 
surement of a distance scale into a constraint on w requires 
various approximations to be made, and depends upon the 
survey in question and upon the time variation assumed for 
the dark energy. Nevertheless it is instructive to go through 
this exercise, bearing these caveats in mind, to get a feel for 
how well future experiments will be able to measure w for 
the case of a constant equation of state. 
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Figure 10. The power spectrum of different galaxy samples measured in real-space, divided by the square of an effective bias parameter 
and the appropriately scaled linear perturbation theory power spectrum. The sample definition and the value of the effective bias used 
are given by the key. The power spectrum of the dark matter spectrum in real-space, also divided by the linear perturbation theory 
spectrum, is shown by the black dashed line. The left hand panel shows the ratios at z = and the right hand panel at z = 1. 
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Figure 11. The ratio of the power spectrum of galaxies measured in redshift-space to that in real-space, at z = (left) and z = 1 
(right). The samples are defined by the key in each panel. The dotted horizontal lines show the predictions of Eq. [5] for the various 
samples. 



The form of the power spectrum of density fluctu- 
ations contains information about basic cosmological pa- 
rameters, and measurements of the galaxy power spectrum 
on large scales have been exploited to ext r act the values 
of th e se parameters ( e .g. ICole et al.l 120051; ISanchez et alj 
20061; iTegmark et ail l200rj ; iPadmanabhan et ail l2007t 
Percival et ail 120071 ). The apparent scale of features in the 
power spectrum offers another route to constrain selected 



cosmological parameters through the dependence of the dis- 
tances parallel and perpendicular to the line of sight on the 
matter density parameter, Qm, the dark energy density pa- 
rameter, Ude, the dark energy equation of state parameter, 
w and the Hubble constant. For such an approach to work, 
we either need to know the true physical scale of a particular 
feature in the power spectrum beforehand or to compare the 
relative size of a feature when measured parallel and perpen- 
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Figure 12. The relation between the dark energy equation of 
state parameter, w, and the scale factor, a, defined by Eq. \E\ 
for perturbations in the equation of state around u>t ruc = 
Two cases are shown. In the upper panel, the values of the other 
cosmological parameters are kept fixed. In the lower panel, the 
ratio of the sound horizon scale to the angular diameter distance 
to the last scattering surface is held fixed. The relation between 
a and w is shown for z = 1 (solid lines) and z = 3 (dashed lines). 
The horizontal and vertical lines guide the eye to show how a 1% 
error in a translates into an error in to. 



dicular to the line of sight (Alcock & Paczynski 1978). The 
baryonic oscillations present a promising candidate for such 
a feature. If we assume for the sake of argument that the 
cosmological parameters, apart from the equation of state 
of the dark energy, are well constrained, then the scale of 
the acoustic oscillations becomes a standard ruler. These 



features are expected on smaller scales than the turnover 
and have already been seen in current surveys at low red- 
shift, although at too low a signal-to-noise ratio to use in 
isolation to extract a competitive constraint on the dark en- 
ergy equation of state (Cole et al. 2005; Eisenstein et al. 
2005). 

We can see how the value of the equation of state param- 
eter parameter of the dark energy influences the form of the 
BAO with the following simple argument. To measure the 
power spectrum of galaxy clustering, we need to convert the 
angular positions and redshifts of the galaxies into comov- 
ing spatial separations. This requires a choice to be made for 
values of the cosmological parameters, including w. In our 
case, we set the parameters equal to the values used in the 
N-body simulations, with w — Wtrue = — 1 for the particular 
case we have run. The effect of a change in the value of w, 
^assumed = TOtruo + Sw is to change the separations between 
pairs of galaxies, which leads to a change in the appearance 
of the power spectrum. For small perturbations away from 
the true equation of state, we assume that the alteration 
in the measured power spectrum can be represented by a 
rescaling of the wavenumber from fc trU c to fc apP . The ratio 
of these wavenumbers gives a "stretch" parameter, a, which 
describes the change in the recovered oscillation scale: 



ktx 



(7) 



If Waaaumed = Wtruej then there is no shift in the BAO in 
the estimated power spectrum and a = 1. In the case of a 
wide-angle, deep galaxy survey with spectroscopic redshifts, 
the stretch parameter can be approximated by: 

n * 1 — ! " 2/3 ( „?(*^™Q , ( 8 ) 



( Da(z, Wassumed) 
Y Da(z, Wtrue) 



H ( Z , Wassumed ; 



where 
H(z,w) 

Da(z, w) 



= Ho [n m (i + zf + n DE (i + z ) 3(1+lu) ] 1/2 

dz 



l + z 



H(z) 



(9) 
(10) 



The values of the exponents in Eq. [8] 2/3 for the distance 
transverse to the line of sight and 1/3 for the distance paral- 
lel to the line of sight are motivated by the number of carte- 
sian components in these directions (e.g. lEisenstein et al.l 
2005). The precise value of these exponents will depend upon 
the geometry and construction of the galaxy survey. For ex- 
ample, in a survey which relies upon photometric redshifts, 
the exponent parallel to the line of sight would be greatly 
reduced and it would be beneficial to compute the power 
spectrum transverse to the line of sight. Note that in Eqs. 9 
and 10 we assume that w is independent of redshift. There 
are many models in which w is a function of redshift. In this 
case, the exponent of £7de in the expression for the Hub- 
ble parameter (Eq. 9) would be replaced by an integral over 
w(z). 

It is instructive to see how the constraints on a translate 
into limits on the value of w. We can do this approximately 
using Eq. [5] for the case of a redshift independent equa- 
tion of state, considering perturbations around u>t ruo = — 1. 
We consider two illustrative cases: a "pessimistic" case in 
which we consider the constraints from BAO in isolation 
from any other data which constrains the cosmological pa- 
rameters and an "optimistic" case, in which we perturb w 
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and only consider cosmological models that give similar pre- 
dictions for the CMB0 The translation in the pessimistic 
case is shown in the upper panel of Fig. [12] for two different 
redshifts. Here we have assumed fixed values for Qm and Qa 
and we have not marginalized over these parameters. This 
is the case discussed most commonly in the literature. Un- 
der these conditions, at z = 1, a 1% error in a corresponds 
approximately to a 4% error in the value of w. At z — 3, the 
boost is about 50% larger, with 8w » 65a. 

In the "optimistic" case, we only consider models which 
give the same angular location for the first peak in the CMB 
spectrum. Hence, when the value of uu is perturbed, we re- 
strict our attention to those models which give the same 
ratio of the sound horizon scale to the angular diameter 
distance to the last scattering surface as our default cosmol- 
ogy. Given the parametric forms quoted for these distances 
by Eisenstein & Hu (1998), this is equivalent to keeping 
Qb/^M and h fixed, and varying Qm- We have called this 
case "optimistic" because it does not include any error on 
the fixed parameters. In this scenario, shown in the lower 
panel of Fig. 1121 the error on w is now only around 50% 
larger than the corresponding error on a. 

We now explore two of the approaches which have been 
advocated in the literature to measure the value of w. Both 
methods involve making fits to the ratio of a measured power 
spectrum divided by a smooth reference spectrum. In the 
first approach, a parametr ic form is assumed for the ratio 
(Blak e fc Glaze brook 20Q3|). The second approach is more 
general as it does not assume a specific form for the ratio, but 
instead uses the linear perturbation th eory power spectrum 
without any further approxim ations (|Percival et al.l 120071 ; 
see also lEisenstein et al.l I2005T ). We shall henceforth refer 
to these methods as the parametric and general schemes re- 
spectively. In their original forms, there are also differences 
in the way in which a "featureless" reference spectrum is 
constructed, as we will briefly discuss when describing these 
approaches below. 

Blake & Glazebrook (2003; see also Glazebrook & Blake 
2005) studied the feasibility of extracting measurements of 
the acoustic oscillations from forthcoming galaxy surveys 
using linear perturbation theory. Their starting point is to 
divide the power spectrum, including the imprint of baryons, 
divided by a smooth reference spectrum which is chosen 
to be free from any signature of acoustic oscillations. This 
method therefore does not use any of the information con- 
tained in the overall shape of the power spectrum, which 
Blake & Glazebrook argue could be susceptible to large 
scale gradients arising from the effects we discussed in Sec- 
tion 3, such as galaxy bias or redshift-space distortions. In- 
stead, they focused on the location and amplitude of the 
acoustic oscillations. The smooth reference spectrum is ob- 
tai ned using the zero-bary on transfer function written down 
by lEisenstein fc Hul (Il99&f ). The parametric form suggested 
by Blake & Glazebrook as a fit to the resulting ratio is a 
Taylor expansion of the ratio of a power spectrum for cold 
dark matter plus a small baryonic component, divided by 
a pure cold dark matter power spectrum. The sound hori- 
zon, which is a free parameter in their method, is treated as 



the oscillation wavelength in this parametric form. This is 
an approximation, as the wavelength of the acoustic oscilla- 
tions actually changes with wavenumber, albeit slowly, and 
is the refore not a constant (see eqn. 22 of lEisenstein fc Hul 
ll998t ). Some authors have criticized this approach due to the 
sensitivity of the ratio to the choice of the reference power 
spectrum. lAngulo et all (|2005l ) describe how realistic power 
spectra, which include nonlinear growth, bias effects and 
redshift-space distortions, require a "linearization" process 
before they become adequately described by the paramet- 
ric form put forward by Blake & Glazebrook. Due to the 
sensitivity of the rat io to the choice of re ference spectrum 
at low wavenumbers, iKoehler et ail (|2006m ) proposed ignor- 
ing power spectrum measurements below fc ~ 0.05MVIpc _1 
to avoid this problem (although we note that they also dis- 
cuss a different approach to measuring the equation of state 
par ameter) . 

iPercival et all (|2007T l proposed a new technique which 
has a number of appealing features compared with that of 
Blake & Glazebrook. Firstly, the shortcut of fitting an ap- 
proximate parametric form to the ratio of the measured 
power spectrum to a reference is dropped in favour of us- 
ing a full linear perturbation theory power spectrum (with 
a modification; see later) to model the ratio. This is com- 
pletely general, and permits one to use the most accurate 
description available of the linear perturbation theory power 
spectrum, such as the tabulated output of CAMB. Secondly, 
the reference power spectrum is defined separately in the 
case of the data and the linear theory model, by using a 
coarse rebinning of the relevant power spectrum. The refer- 
ence is constructed using a spline fit to a reduced number 
of wavenumber bins over the range in which the spectrum 
in question is defined. Thus, any deviations in the general 
form of the measured spectrum away from linear theory are 
naturally accounted for in the reference spectrum. Thirdly, 
Percival et al. allow for a damping of the amplitude of the 
oscillations in the theoretical ratio beyond some wavenum- 
ber, which is treated as a free parameter in their fit. The 
quality of the fits is dramatically improved when damping 
of the higher harmonics is allowed. Percival et al. applied 
their method to extract the matter density parameter from 
the power spectrum of luminous red galaxies in the SDSS. 

The majority of the results we present are obtained us- 
ing the general method suggested by Percival et al. For com- 
pleteness, and because Percival et al. did not actually apply 
their method to the extraction of the equation of state pa- 
rameter, we set out the general approach step-by-step below: 

1. A smooth reference spectrum (i.e. without any oscilla- 
tory features), P rc f, is constructed from the measured power 
spectrum using a cubic spline fit over the wavenumber range 
0.0046 < (fc/ZiMpc -1 ) < 1.2, using the measured spectrum 
smoothed over 25 bins in wavenumber. The spline is con- 
strained to pass through the data points in this coarse re- 
binning of the measured power spectrum. 

2. We compute the ratio, R(k), of the measured power 
spectrum, P(k), to the reference spectrum, P rG f (fc), obtained 
in step 1: 



2 We acknowledge the referee for suggesting this second case to 
us and for encouraging us to perform the calculation. 



(11) 



3. A linear perturbation theory power spectrum is gener- 
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ated with CAMB for the cosmological parameters used in the 
BASICC simulation. A smooth reference spectrum, P^f j is de- 
fined for this spectrum in the same manner as described for 
the measured spectrum in Step 1, using the same wavenum- 
ber bins. A ratio, J?l, is derived for the linear perturbation 
theory spectrum by dividing by this reference spectrum. 

4. The linear theory ratio, _Rl, is compared with the mea- 
sured ratio, R. Two modifications are considered to the 
linear theory ratio. The first is a stretch or scaling of the 
wavenumber used in the linear theory ratio, as described 
above, to mimic the act of changing the dark energy equa- 
tion of state parameter, w. The goal here is to see what 
variation in w can be tolerated before -Rl is no longer a 
good fit to the measured ratio R. The second change is to 
allow for a damping of the oscillations beyond some char- 
acteristic wavenumber by multiplying the theoretical power 
spectrum by a Gaussian filter: 



W(k) = exp 



(12) 



where k n i is a free parameter. Hence, the linear theory ratio 
is modified to: 



(13) 



RUk)= ( pr(o*)-l ) xW(k,k nl ) + l 



5. A likelihood is computed for each combination of the 
parameters fe n i and a, assuming Gaussian errors: 



21nL = x = 



R'-Rj 
a i jP i 



(14) 



where the summation is over wavenumber and a 1 is the error 
on the power spectrum estimated in the i th bin (as given by 
Eq. [3}. We generate a grid of models using 200 2 different 
combinations of a and k n i in the ranges [0.9,1.1] and [0,0.4] 
respectively. 

6. Finally, the best fit values for a and k n \ correspond to 
those for the model with the maximum likehood. We obtain 
confidence limits on the parameter estimation by considering 
the models within A^ 2 equal to 2.3 and 6.0; in the case of 
a Gaussian likelihood, these would correspond to the 68% 
(l-cr error) and 95% (2-a error) confidence levels on the best 
fit. We note that in some cases presented later (see Fig. 15), 
the distribution of the likelihood is not Gaussian. 

In some cases, we also present constraints on w derived us- 
ing a slightly modified version of the approach of Blake & 
Glazebrook. The main difference is that we follow step 1 to 
construct a ratio from the measured power spectrum, rather 
than using a zero-baryon transfer function. 

One issue to be resolved is the range of wavenumbers 
which should be used in the fitting process. To address this, 
we used the power spectrum of the dark matter measured 
at z — 6. We systematically varied the minimum and max- 
imum wavenumbers used in our fit and compared the val- 
ues of the scaling parameter, a, recovered. Our results are 
fairly insensitive to the choice of the maximum wavenumber, 
particularly when damping of the oscillations is included 
in the fitting algorithm. However, the recovered a shows 
a systematic shift once the minimum wavenumber exceeds 
k ~ 0.1/iMpc -1 . For minimum wavenumbers smaller than 
this, there is little difference in the recovered value of a 
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Figure 13. The power spectra of dark matter particles, dark 
matter haloes and galaxies at z = (error bars). The real- 
space power spectra are plotted in the left hand column and the 
redshift-space power spectra appear in the right hand column. 
The red curves show the reference spectra derived from the mea- 
sured spectra using a cubic spline fit, as described in Section 4. 
The blue curve is the same in each panel, showing the linear per- 
turbation theory prediction for the 2 = matter power spectrum 
(plotted using a high redshift output obtained from the BASICC 
simulation, which has been scaled in amplitude according to the 
difference in growth factors between the two epochs expected in 
linear perturbation theory) The errors on the power spectrum are 
estimated using Eq. [3] 



or in the size of the errors on a, as these modes have rel- 
atively large errors in our simulation. This is encouraging 
news for realistic survey geometries, for which the power 
spectrum measured at low wavenumbers will be distorted 
due to the window function of the survey. In the rest of 
the paper, we use the power spectrum in the wavenumber 
interval fc/(MVIpc -1 ) = [0, 0.4] to constrain the value of a. 



5 RESULTS 

In this Section, we present the expected constraints on the 
dark matter equation of state using the power spectra mea- 
sured from our simulations. We first show how our algo- 
rithm for extracting the equation of state parameter works 
in practice, for dark matter particles, haloes and galaxies, 
comparing the results obtained in real-space and redshift- 
space (§5.1). We then assess the need for an accurate model 
of the linear theory power spectrum and the relative merits 
of the general and parametric fitting procedures (§5.2). In 
§5.3, we present our main results, which are summarized in 
Fig. [19] and Table [2] which lists the best-fitting value of a 
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Figure 14. The ratio of the measured power spectrum divided 
by a smooth reference spectrum. The symbols correspond to the 
measurements plotted in Fig. [13] divided by the red curve in each 
panel of that figure. The red lines here show the best- fitting model 
in each case using the general method and the blue curves show 
the best fit for the parametric method. The errors on the power 
spectrum are estimated using Eq. [3] 



Figure 15. The constraints on the parameters fc n i and a for the 
power spectra plotted in Fig. 1131 The contours show the 1, 2 and 
3-(T confidence limits for two parameters. 



and the estimated error for different samples of galaxies at 
z = 1, along with the corresponding fractional error in w. 
Finally, in §5.4, we use the results presented in §5.3 to make 
forecasts for the accuracy with which several forthcoming 
surveys will be able to measure the value of w. 



5.1 The algorithm to extract the scale of the 
acoustic oscillations in action 

We present a series of plots for samples at z = 0, which 
illustrate the various stages in the fitting process. Fig. 1131 
shows the power spectra measured for different tracers, both 
in real-space and redshift-space. The sample of dark mat- 
ter haloes includes all objects with a mass in excess of 
5.4 x 10 12 Mq. The galaxy sample is magnitude-limited 
with a space density of ft = 5 x 10~ 4 h~ 3 Mpc 3 . For ref- 
erence, the linear perturbation theory power spectrum for 
the mass at z = is shown by the blue line in each panel: 
this is the power spectrum of the dark matter measured in 
real-space at z = 15, scaled by the ratio of growth factors 
in order to have the amplitude expected at z = 0. It is 
important to bear in mind that the y-axis in this plot cov- 
ers more than a factor of one thousand in amplitude. Fig. 1131 
shows that there is considerable variation in the power spec- 
tra measured for different types of objects, and between the 
results in real-space and redshift-space, which re-inforces the 



points made in Section 3 regarding deviations from the pre- 
dictions of linear perturbation theory on large scales. The 
red curve in each panel shows the corresponding reference 
power spectrum, which is constructed from the measured 
power spectrum as explained in Section 4. 

In Fig. 1141 the symbols show the ratio obtained by di- 
viding the measured power spectrum by the appropriate ref- 
erence spectrum for the same samples plotted in Fig. 1131 
The ratios look remarkably similar for the different trac- 
ers up to k ps 0.15/iMpc -1 . Beyond this wavenumber, the 
appearance of the oscillations varies from panel to panel, 
but the ratio stays close to unity. This similarity illustrates 
how well the approach for producing the reference spectrum 
works. The red curves in each panel show the best-fitting 
model produced in the general scheme whilst the blue curves 
show the fit obtained in the parametric approach. The best 
fits have somewhat different forms at wavenumbers below 
k ~ 0.05ftMpc _1 . The constraints on the values of the pa- 
rameters fc n i and a are presented in Fig. 1151 where we show 
the 1, 2 and 3-cr ranges in the case of two parameters, com- 
puted assuming Gaussian errors. There is a weak systematic 
trend for the best-fitting result for a to shift to lower val- 
ues when galaxies are considered instead of the dark matter. 
The errors on the recovered parameters are larger in the case 
of galaxies than for the dark matter or for haloes, reflect- 
ing the lower signal-to-noise of the predicted galaxy power 
spectrum. 
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Figure 16. The best fit value for the scaling parameter a, re- 
covered from the ensemble of low resolution simulations, using 
the dark matter power spectrum in real-space. The results are 
show for two different redshifts: z = 3.8 (top) and z = bottom. 
The histograms marked CAMB and BG03 show the results for the 
general and parametric fitting procedures, respectively. The blue 
histogram shows the results if the general method is followed with 
the CAMB power spectrum replaced by the formula for the linear 
theory power spectrum presented by Eisenstein & Hu (1998). 



5.2 Two tests of the algorithm 

Before presenting the main results of applying our algorithm 
to extract the acoustic oscillation scales for various samples 
drawn from the BASICC run, we use the L-BASICC ensem- 
ble to address two questions: 1) How accurately do we need 
to model the linear perturbation theory matter power spec- 
trum to avoid introducing a systematic bias into the results 
for the oscillation scale? 2) How does the performance of the 
new method for constraining the oscillation scale introduced 
in this paper compare with earlier approaches? To help an- 
swer these questions, we use the power spectrum of the dark 
matter measured from the L-BASICC runs in real-space at 
z = and z = 3.8, the highest output redshift besides the 
initial conditions. The results of applying our standard al- 
gorithm for extracting the oscillation scale are shown by the 
red histogram labelled CAMB in Fig. 1161 which gives the dis- 
tribution of the best-fitting value of a. The ensemble returns 
an unbiased mean value for the stretch parameter, a = 1. 



At z = 3.8, the standard deviation on the best fit is 0.3%; 
by z = 0, this rises to 1%. 

To address the first issue above, regarding how well we 
need to model the linear theory power spectrum to get an 
unbiased result for the oscillation scale, we replace the CAMB 
generated power spectr um in our algorit h m by the approx- 
imation introduced by lEisenstein fc Hul (|l998l ). These au- 
thors proposed a physically motivated expression for the lin- 
ear theory power spectrum, with parameters set to achieve a 
reasonable match to the results obtained from detailed cal- 
culations using Boltzmann codes over a much wider range 
of wavenumbers than are typically considered for baryonic 
acoustic oscillations. Eisenstein & Hu's motivation was to 
provide physical insight into the form of the power spec- 
trum in a cold dark matter universe and to produce a code 
which could rapidly calculate large numbers of power spec- 
tra for grids cosmological parameters. Of course, the correct 
approach in our fitting procedure is to use the same code to 
compute the linear theory spectrum as was used to gener- 
ate the initial conditions in the N-body simulation. In the 
case of real data, we do not have the luxury of knowing 
which Boltzmann code to use, so we should use the one 
which claims to be the most accurate representation of the 
model we are testing. Nevertheless, it is instructive to per- 
form this test to see what error is introduced by using a less 
accurate calculation of the transfer function. The choice of 
Eisenstein & Hu's code is particularly relevant for this pur- 
pose as Blake & Glazebrook used this formalism to inspire 
their parametric expression to fit the acoustic oscillations. 
The use of Eisenstein & Hu's formalism to model the lin- 
ear theory power spectra generated with CAMB introduces a 
small but measurable systematic shift in the mean value of 
a. At z = 0, the mean a indicated by the blue histogram in 
Fig. [16] is 0.98 ±0.01. 

We answer the second question by adopting the fitting 
algorithm of Blake & Glazebrook (2003), which assumes a 
parametric form for the ratio of the power spectrum with 
baryons to a smooth, cold dark matter only power spectrum. 
Changing the fitting method in this way also introduces a 
similar magnitude of shift in the best-fitting value of a. The 
green histogram shows the results when we use the para- 
metric approach introduced by Blake & Glazebrook (2003). 
The mean value of a in this case is 1.01 ± 0.01. These shifts 
are small but one must bear in mind that the correspond- 
ing bias in the dark energy equation of state parameter is 
several times larger than the shift in a. 

dsf 



5.3 The main results 

We now turn our attention back to the general results shown 
in Fig. [17] and Fig. 1181 and discuss the conclusions for dif- 
ferent tracers of the density field in turn. In these plots, 
the symbols refer to the constraints obtained from the high 
resolution BASICC simulation and the shading shows results 
from the ensemble of low resolution simulations, L-BASICC. 

The blue triangles in Fig. [T7] show the values obtained 
for a from the power spectrum of the dark matter. There 
is a trend for the best-fitting value to deviate away from 
unity with decreasing redshift, although the result at z — 
is still within 1-a of a — 1. The mean of the ensemble of low 
resolution runs does not, however, show any deviation away 
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Table 2. The results of applying the general fitting procedure described in §4 to power spectra measured for different galaxy catalogues 
at z = (top) and z = 1 (bottom). In each table, the first row gives the results for the dark matter and the final row lists results for a 
sample of dark matter haloes (all haloes with mass in excess of 2.7 X 10 13 h' 1 M Q ). The first column gives the label of the sample, as 
defined in Section 2. The second column gives the space density of galaxies. The first two samples, A and B, are constructed by applying a 
magnitude limit. Samples C-F are derived from sample A by applying a second selection criterion, as listed in the third column. Samples 
C and E correspond to the red and blue halves of sample A respectively. Samples D and F comprise the 50% of galaxies from sample 
A with the strongest and weakest (in terms of equivalent width) OII[3727] emission lines, respectively. Column 4(10) gives the effective 
bias of the sample, computed from the square root of the ratio of the measured galaxy power spectrum in real (redshift) space to the real 
space power spectrum of the dark matter over the wavenumber interval 0.01 < (fc//iMpc -1 ) < 0.05. Column 5(1') gives the ratio of the 
clustering signal to the shot noise for the power spectrum measurement, averaged over the wavenumber range 0.19 < (fc//iMpc _1 ) < 0.21. 
Columns 6 and 7 (12 and 13) give the best-fitting values of the scaling parameter a and the 1 — cr error on the fit, in real (redshift) space. 
Column 9 (15) gives the error expected on the scale parameter from Seo & Eisenstein (2007). The rms lagrangian displacement was set 
equal to 1 over the best fit non linear scale (l/fc nI ) for each case. 



from a = 1 as a function of redshift, although the scatter on 
the recovered value of a increases towards the present day. If 
we examine the analogous results for individual simulations 
taken from the low resolution ensemble, we find a wide range 
of behaviour for the best-fitting value of a for the dark mat- 
ter. Some low resolution runs give results which look like the 
high resolution one, whereas others show deviations away 
from a — 1, with values of a > 1, as z = is approached. 
The trend seen for the dark matter in the high resolution run 
serves to illustrate the importance of sampling fluctuations, 
even in such large volumes. In redshift-space, the scatter in 
the recovered value of a is larger than in real- space (see also 
ISeo fc Eisenstein! 120051 : lEisenstein et al.ll2006J ). 

To obtain the errors quoted in Table 2 on the param- 
eters a and k n \, we assume Gaussian mode counting errors 
on the power spectra measured in the BASICC simulation, as 
given by Eq. 3. In Fig. 3, we showed that this simple estimate 
of the errors on the power spectrum agreed fairly well with 



the scatter found in the measurements from the L-BASICC 
ensemble, particularly for the case of the dark matter. We 
have extended this comparison to look at how the errors 
on a and fc n i quoted in Table 2 match the scatter in these 
parameters obtained from the L-BASICC runs. We find the 
scatter estimated from the ensemble is somewhat larger than 
the error inferred using the mode counting argument. At 
z = 0, the mode counting errors are 20% smaller for a for 
the dark matter in real space. In redshift-space, the discrep- 
ancy increases to nearly 30%. The mismatch between the 
two estimates is smaller at z = 1. The level of disagreement 
is not remarkable. It could be the case that the scatter from 
the ensemble has not converged, even with 50 realizations of 
the density field. A more likely explanation, particularly in 
view of the redshift dependence of the discrepancy, is mode 
coupling in the power spectrum measurements arising from 
nonlinearities and redshift space distortions, which could in- 
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Figure 17. The best-fitting value of the scale factor a as a 
function of redshift, for different tracers of the density distribu- 
tion, in real-space (top) and redshift-space (bottom). The sym- 
bols show results from the high resolution BASICC simulation: dark 
matter (blue triangles), dark matter haloes with mass in excess 
of 5.4 X 10 12 h' 1 M e ( green circles) and galaxies (red squares). 
The error bars show the 1-cr range on a, calculated from Ax 2 . 
The hatched region shows the central 68% range of the results 
obtained using the dark matter in the ensemble of low resolu- 
tion simulations. Recall that a = 1 corresponds to an unbiased 
measurement of the equation of state parameter, w, and that 
8w « 45a at z = 1. 



crease the variance in the power spectrum compared with 
the Gaussian estimate. 

Fig. [18] shows that there is a strong trend for the best- 
fitting value of the smoothing scale, fc n i, to decrease with 
decreasing redshift. This results from the oscillations being 
erased and modified down to smaller wavenumbers as the 
nonlinearities in the density field grow. The variation of the 
smoothing scale k n \ on redshift is well described by a linear 
relation: fe n j — a + bz. In real-space, a = 0.108 ± 0.0082 and 
b = 0.054 ±0.01 10. In redshift-space, a = 0.096 ±0.0074 and 
b = 0.036 ± 0.0094. 

The constraints on a and fc n i for dark matter haloes 
(with masses in excess of 5 x 10 12 h -1 Mq) axe plotted with 



Figure 18. The best-fitting value of the damping scale 

"'111 ^ 

function of redshift, for different tracers of the density distribu- 
tion, in real-space (top) and redshift-space (bottom). The sym- 
bols show results from the high resolution BASICC simulation: dark 
matter (blue triangles), dark matter haloes with mass in excess of 
5.4 X 10 12 h" 1 Mq (green circles) and galaxies (red squares). The 
error bars show the 1-<t range on a. The hatched region shows the 
central 68% range of the results obtained using the dark matter 
in the ensemble of low resolution simulations. 



green circles in Figs. [T7] and 1181 The parameter constraints 
obtained for this sample of haloes are very similar to those 
found for the dark matter, except for the value of k n \ at 
high redshift. Considering haloes in place of dark matter 
represents a step closer to the observations, so it is reassuring 
that the conclusions do not change significantly. 

Finally, in Figs. [17] and [18] we show using red squares 
the results for magnitude-limited samples of galaxies. The 
magnitude limit is varied with redshift such that in each 
case the galaxy sample has a space density n = 5 x 
10 -4 h~ 3 Mpc 3 . There is a weak systematic shift in the best- 
fitting values of a compared with the results obtained for 
the dark matter. At the same time the signal to noise of the 
power spectrum measurement is lower for the galaxy samples 
than for the dark matter, so the errors on the best-fitting 
parameters are correspondingly larger for the galaxies. The 
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Figure 19. The recovered value of the stretch parameter a for 
the galaxy samples listed in Table 2. Recall that a = 1 corre- 
sponds to the equation of state parameter w = — 1. At z = 1, a 
shift in a away from unity implies a shift in the recovered value 
of w given by Sw ~ 45a. 



galaxy samples are consistent with a — 1 at slightly over 1-a. 
The size of this systematic shift is comparable to the random 
measurement errors, so we cannot reach a firm conclusion. 
It will be very interesting to repeat our calculation with a 
larger simulation volume to reduce the size of the random 
errors and to assess if such shifts could genuinely provide an 
ultimate limitation to the accuracy of this method. 

As a result of using a semi-analytic galaxy formation 
model which makes predictions for the observable proper- 
ties of galaxies, we can vary the selection criteria used to 
construct samples and compare the constraints on the equa- 
tion of state. The results of this exercise at z = 1 are pre- 
sented in Table 2 and in Fig. 1191 where we consider a range 
of samples defined either by a simple magnitude limit, or 
by a magnitude limit applied in combination with a colour 
cut or a restriction on the strength of an emission line. The 
key result from comparing the constraints for different sam- 
ples is that whilst there are no strong systematic differences 
between the results, the accuracy of the constraints varies 
significantly. For example, using a catalogue of red galaxies, 
we predict that one could measure the dark energy equation 
of state with an accuracy 40% better than with the same 
number density of galaxies chosen by the strength of their 
emission lines. 

We compare the error on the acoustic scale extracted 
from our s imulations with the res ults of the prescription 
set ou t by ISeo fc Eisensteinl (|2007T) . The ISeo fc Eisensteinl 
(|2007l ) algorithm contains a parameter which is equivalent 
to 1/fcni- If we use our best fitting values of fc n i, we find that 
the Seo & Eisenstein prescription gives a similar estimate of 
the error on the acoustic scale to that we obtain by fitting 
directly to the simulation results. However, if we use the 
value of fc n i suggested by ISeo fc Eisensteinl ()2007h ■ which 



they extract from a dark matter simulation, we find that 
their prescription gives an optimistic estimate of the error 
on a. The reason we recover a larger value of k n \ from our 
galaxy samples than we do for the dark matter is due to 
the increased discreteness shot noise in these samples, which 
results in noisier power spectra at high k. This causes an 
elongation in the confidence levels in the kni versus a plane. 

It is interesting to compare the results for the dark mat- 
ter and for the galaxy samples with those for a set of massive 
haloes. Table 2 also gives the constraints on a and k n i for a 
sample of massive haloes (see also Angulo et al. 2005) . There 
are 142 000 haloes in the BASICC output at z=l with a mass 
in excess of 2.7 x 10 13 h~ x Mq. Although the effective bias of 
this sample of massive haloes is greater than that of any of 
the galaxy samples listed in Table 2, the reduction in space 
density means that nP ~ 1 and the estimated error on w is 
comparable to that found for the galaxy samples. 

5.4 Forecasts for future surveys 

We can use the results presented in Table 2 to make a rough 
estimate of the accuracy with which future surveys are likely 
to be able to constrain the scale of the acoustic oscillations. 
This can be done using a simple calculation motivated by 
the expression for the fractional error in the power spectrum 
given by Eq. 3. We assume that the error in the distance 
scale, Aq, scales with the volume of the survey, V^ urV oy, and 
the product of the space density of galaxies and the power 
spectrum, nP(k = 0.2/iMpc - ), as: 

1 



Aq ex 



(15) 



The constant of proportionality can be set for a particular 
galaxy sample using the results given in Table 2. 

The WiggleZ survey is currently underway and will 
measure redshifts for 400,000 blue galaxies over 1000 square 
degrees in the redshift interval z = 0.5 — 1.0 (Glazebrook 
et al. 2007). For the cosmological parameters adopted in 
this paper, this gives a comoving volume of 1.13 h~ 3 Gpc 3 . 
Using the blue colour selected sample or the large equiva- 
lent width sample from Table 2, and assuming nP ~ 1 for 
WiggleZ galaxies, somewhat higher than we find in our sim- 
ulation, we estimate that this survey will measure the dis- 
tance scale to an accuracy of Acv ~ 2%, which is similar to 
that claimed by Glazebrook et al. using linear perturbation 
theory. 

The WFMOS survey has been proposed to motivate the 
construction of a new spectrograph for the Subaru telescope 
(Glazebrook et al. 2005). This will target galaxies with a 
space density of n — 5 x 10 -4 h 3 Mpc -3 in the redshift in- 
terval z — 0.5 — 1.3 over 2000 square degrees, covering a 
volume of 4.4 /i -3 Gpc 3 . (There is also a WFMOS survey 
which will target z — 3 galaxies but over a much smaller 
solid angle.) Using sample A from Table 2, and adopting 
nP — 1, we obtain an estimated error of Aa = 0.83%, again 
in good agreement with Glazebrook et al. 

Photometric surveys can generally cover a larger solid 
angle than spectroscopic surveys down to a fainter magni- 
tude limit. The fainter magnitude limit results in a higher 
median redshift and a broader redshift distribution for the 
survey galaxies, which means that a larger volume is cov- 
ered. However, the limited accuracy of photometric redshift 
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estimates means that in practice Fourier modes are lost and 
the effective volume of the survey is greatly reduced. Blake 
& Bridle (2005) estimate that the factor by which the survey 
volume is reduced is « 12 (Sz/(1 + z)/0.03), where 5z/(l + z) 
is the error in the photometric redshifts. 

The Panoramic Survey Telescope and Rapid Response 
System (Pan-STARRS) survey will map 3n steradians of the 
sky ( |http: / / pan-st arrs.ifa.hawaii.edu /public /h ome.html[ ) . 
Cai et al. (2007, in preparation) show that the median 
redshift of the 3tt survey will be z ~ 0.5, with a tail 
extending to z ~ 1.2. The volume of the survey, assuming 
that 20, 000 square degrees cover low-extinction parts of 
the sky and give high quality clustering measurements, is 
around 41/i _3 Gpc 3 . Talking sample A from Table 2, and 
setting nP >> 1, as appropriate for the relatively high 
space density of galaxies in a photometric sample, and 
allowing for the reduction in the effective volume caused by 
a photometric redshift error of Sz/(1 + z) = 0.03, gives a 
forecast error on the oscillation scale of Aa ~ 0.5%. In the 
more likely event that the photometric redshift errors are 
twice as large, Sz/(l + z) ~ 0.06, this figure increases to 
Aa ~ 0.7%. 

Remembering the crude conversion Aw » 4Aa from 
Section 4, this means that the next generation of galaxy 
surveys is unlikely to deliver 1% errors on a constant equa- 
tion of state from BAO measurements used in isolation from 
other cosmological data. A survey with almost an order of 
magnitude more effective volume than Pan-STARRS will be 
needed to achieve this target. This will require an all-sky, 
spectroscopic galaxy redshift survey, such as the SPACE mis- 
sion being proposed to ESA's Cosmic Vision call. SPACE will 
measure redshifts for galaxies in the interval 0.5 < z < 2, 
covering around 150/i _3 Gpc 3 . Extrapolating from Sample 
A, we forecast that an error in the oscillation scale of 
Aa ~ 0.15% could be achieved with SPACE. In the case of 
the pessimistic translation to an error on w considered in 
Section 4, this corresponds to Aw ~ 0.6%; in the optimistic 
scenario, we expect a constraint of Aw ~ 0.23%. 



6 CONCLUSIONS 

In the next five to ten years, several proposed galaxy surveys 
will allow high precision measurements of the clustering of 
galaxies on the scale of the acoustic oscillations at intermedi- 
ate and high redshifts. Both photometric and spectroscopic 
surveys are planned, which will cover volumes up to tens 
of cubic gigaparsecs and will contain hundreds of thousands 
to hundreds of millions of galaxies. There is a clear need to 
ensure that theoretical predictions develop apace with suf- 
ficient accuracy and realism to allow such datasets to be 
fully exploited and to uncover any possible systematic er- 
rors in this cosmological test to uncover the nature of the 
dark energy. 

Early theor etical work in this area use d linear pertur- 
bation theory llBlake fc Glazebrookl 120031 ; iHu fc Haimanl 
120031 : iGlazebrook fc Blake! l2005h . Recently, more physical 
calculations have been carried out using N -body simulations 
with cubes of side 500 - 1100 h ' 1 Mp c dSeo fc Eiscnstein 



20031. 120051: ISchulz fc White! 120061 : iHuff et al l l200d 



Seo fc Eisensteir] 120071 ) . In this paper, we have improved 



upon previous modelling work in three ways. Firstly, we 



have used a simulation volume comparable to the largest of 
the currently proposed spectroscopic surveys. This allows 
us to accurately follow the growth of density fluctuations 
on an ultra large scales in excess 100 Mpc, the scales 
of interest for the acoustic oscillations, which can only be 
followed approximately in smaller computational volumes. 
In particular, a large volume is necessary to obtain accurate 
predictions for bulk flows, which are sensitive to the power 
spectrum at low wavenumbers. The only published work 
with a larg er simulation volume used the Hubble V olume 
simulation (| Angulo et all 120051 ; iKoehler et all l2006h . The 
Hubble Volume has a larger particle mass than the BASICC, 
which restricted these studies to consider either cluster mass 
dark matter haloes (Angul o et al. 2005) or a sim ple biasing 
scheme to add galaxies (jKoehler et al] l2006h . Secondly, 
through the use of a large number of particles, we are able 
to resolve the majority of the haloes which are likely to 
host the galaxies which will be observed in the forthcoming 
surveys. Thirdly, we use a semi-analytic galaxy formation 
model to populate the simulation with galaxies. Unlike 
other studies which use phenomenological biasing schemes 
or the halo occupation model to add galaxies, this allows 
us to predict the shape and amplitude of the galaxy power 
spectrum and the signal-to-noise of the clustering expected 
for different galaxy selections. 

We use our N-body simulation in combination with a 
galaxy formation model to make the connection between the 
linear perturbation theory prediction for the matter power 
spectrum and the power spectrum of galaxies. We do this in 
a series of steps, starting with power spectrum of the dark 
matter, looking at the impact of the nonlinear growth of 
fluctuations and peculiar motions or redshift-space distor- 
tions, before examining the power spectrum of dark mat- 
ter haloes and, finally, galaxies. A number of conclusions 
are reached from this study: i) The nonlinear evolution of 
the dark matter power spectrum is apparent even on scales 
larger than the sound horizon scale. Although the devia- 
tion from linear theory is only a few percent, the coupled 
evolution of different Fourier modes means that these scales 
need to be followed accurately to get the correct behaviour 
at higher wavenumbers. ii) The form of the distortion of 
the power spectrum due to peculiar motions is extremely 
sensitive to the type of object under consideration, being 
quite different for the cases of dark matter, dark haloes and 
galaxies. Moreover, different galaxy selections give different 
redshift-space distortions, iii) Galaxy bias is scale depen- 
dent and sensitive to the selection applied for wavenumbers 
k > 0.15ftMpc~ . Eisenstein et al. (2006) discuss a technique 
which attempts to reconstruct the linear density field from 
an observed distribution of objects. The reconstruction can 
reduce the damping of the higher harmonic oscillations in 
the power spectrum, thereby increasing the statistical sig- 
nificance of the acoustic scale measurement and diminishing 
any systematic effects caused by departures from linearity. 
It will be interesting to apply this method to the galaxy 
samples presented in this paper, to see if this approach still 
works at the required level in the case of biased tracers of 
the linear density field. 

We also present a new method to extract the dark en- 
ergy equation of state parame t er, ba sed upon an approach 
put forward bv lPercival et all (|2007T ) . The method involves 
dividing the measured power spectrum by a smooth refer- 
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ence spectrum and comparing the resulting ratio to the pre- 
dictions of linear perturbation theory. The algorithm has 
th ree key advances over e arlier work, which can be credited 
to lEisenstein et al.1 (120051 ) and lPercival et al.1 (|2007h : i) The 
reference spectrum is derived from the measured spectrum, 
which avoids the need to apply major corrections to a lin- 
ear theory reference, ii) The measured ratio is compared to 
a prediction generated using CAMB, which is more accurate 
than assuming a parametric form for the ratio based on a 
Taylor expansion, iii) The linear theory ratio is modified by 
allowing the higher-order oscillations to be damped, which 
improves the fit to the measured ratio. Changing the value of 
the equation of state parameter is approximately equivalent 
to rescaling the wavenumber in the predicted power spec- 
trum ratio; at z = 1, a 1% shift in wavenumber is equivalent 
to a 4% shift in the recovered value of w. 

We explore the constraints on the dark energy equation 
of state using different tracers of the density field. By ap- 
plying our algorithm for extracting the oscillation scale to 
the L-BASICC ensemble, we have provided the most stringent 
test to date of usefulness of baryonic acoustic oscillations for 
measuring the equation of state of the dark energy. For the 
case of the dark matter, there is no significant bias in the re- 
covered oscillation scale, compared with the value expected 
from linear perturbation theory. Within a given simulation, 
we find that 1% deviations from the underlying length scale 
are possible although these are only at the 1-a level. Such 
excursions are the result of sampling variance arising from 
the finite volume of the computational box, which are im- 
portant even in a simulation of the volume of the BASICC. 
The error on the scale factor recovered from galaxy samples 
is larger than that found for the dark matter, reflecting the 
lower signal-to-noise of the galaxy power spectrum measure- 
ments. Different galaxy selections lead to variations in the 
clustering strength and hence in the error expected in the 
scale factor. 

Currently, the best constraints on the equation of state 
parameter come either from combining datasets, such as 
the power spectrum of galaxy clustering and measurements 
of the microwave background radiation (e.g. Sanchez et al. 
2006) or from the Hubble diagram of Type la, with priors 
on the flatness of the Universe and the matter density (Riess 
et al. 2004). For example, Wood-Vasey et al. (2007) combine 
high redshift SNe la from the ESSENCE Supernova Survey 
with the measurement of the BAO made by Eisenstein et al. 
(2005), and, assuming a flat universe, constrain a constant 
equation of state to have w = — 1.05 jlj}'^ (stat.)±0.11(sys.), 
consistent with a cosmological constant. Possible contribu- 
tions to the systematic error include the degree of dust ex- 
tinction in the SNe host galaxy, evolution in the properties 
of SNe with redshift and local calibration effects such as a 
"Hubble bubble". We have used our simulation results to 
forecast the accuracy with which future galaxy surveys will 
use the BAO in isolation to constrain the scale of the acous- 
tic oscillations, and under certain assumptions, w. We an- 
ticipate that Pan-STARRS, with accurate photometric red- 
shifts, will have an accuracy comparable to that expected 
for the next generation of spectroscopic survey (WFMOS) 
and could potentially reduce the statistical errors on the 
value of to by a factor of 2 compared with the current con- 
straints. However, the target of 1% random errors on w using 
BAO measurements is beyond the grasp of any of the sur- 



veys likely to be completed or even to start within the next 
decade. 

The predictions we have presented here are idealized 
in a number of respects. The accuracy with which we ex- 
pect the dark energy equation of state parameter will be 
measured assumes that the values of the other cosmologi- 
cal parameters are known with infinite accuracy. We have 
also neglected the impact of the survey window function on 
the power spectrum measurement; this will be particularly 
important in the case of surveys which rely on photometric 
redshifts. In future work, we plan an number of improve- 
ments: i) Use of an even larger simulation volume, to exceed 
that proposed in forthcoming surveys. One caveat on our 
quoted error on w is that some of the planned surveys will be 
larger than the volume of the BASICC, and will consequently 
have smaller sampling fluctuations, ii) The inclusion of the 
evolution of clustering along the line of sight. Although we 
have focused on z = 1, proposed surveys will span a broad 
redshift interval centred on this value, iii) The inclusion of 
a survey window function, mimicing the angular and radial 
selections, and including the impact of errors on photomet- 
ric redshifts. Such calculations represent huge challenges in 
computational cosmology, due to the volume coverage and 
mass resolution required in the N-body simulations used, 
and the post-processing needed to include galaxies. How- 
ever, such calculations are essential if the BAO approach is 
to be used to its full potential. 
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